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This  paper  brings  together  different  aspects 
of  the  Critical  Path  Analysis  exclusively  in  terms 
of  the  mathematical  developments  in  this  area  during 
the  last  eight  years.  The  survey  is  divided  into 
seven  sections.  The  first  one  deals  with  the  four 
stages  of  development,  each  successive  stage  being 
more  representative  of  the  real  situation.  The  next 
five  sections  give  explicit  reference  to  the  different 
articles  which  contribute  to  the  stages  of  development 
described  in  the  first  section.  The  articles  are  re¬ 
viewed  in  terms  of  the  network  algebra,  PERT  assump¬ 
tions,  cost/time  trade  off,  resource  allocation 
(heuristic)  and  resource  allocation  (analytical). 

Finally,  the  seventh  and  last  section  of  this  paper 
provides  a  summary  of  ideas  and  possible  directions 
of  research,  some  of  which  are  currently  being  undertaken. 
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SECTION  I:  STAGES  OF  DEVELOPMENT 


1.1  FIRST  STAGE 

The  origin  of  the  network  approach  to  project 
management  can  be  trac:  i  back  to  two  separate  projects; 
one  undertaken  by  the  private  sector  and  the  other  by 
the  U.  S.  Government,  the  latter  being  very  much  larger 
in  magnitude  and  in  the  degree  of  uncertainty  than  the 
others.  Both  groups  advocated  the  use  of  a  network  de¬ 
picting  explicitly  the  relationship  among  various  ac¬ 
tivities.  This  was  indeed  a  significant  change  from 
the  then  existing  practice  of  using  barcharts.  The 
smaller  of  the  two  projects  related  to  the  construction 
industry  where,  by  virtue  of  extensive  experience,  the 
activities  could  be  well  defined  at  the  planning  stage 
and  unique  time  values  ascribed  to  each  one.  This  type 
of  network  can  be  classified  as  deterministic  and  activ¬ 
ity  oriented.  Kelley  and  Walker  [64]  had  initiated  this 
under  the  title  "The  Critical  Path  Method". 

The  method  consists  of  drawing  a  network  of  the 
project  portraying  the  interrelationship  among  the  ac¬ 
tivities  and  assigning  time  values  (dij)  to  each  activ¬ 
ity  { t,j  ) .  Every  activity  is  uniquely  represented  by 
a  predessor  event('-)  and  a  successor  eveni(j)  such 
that  L<J. 


4. 


The  time  calculations  involve  a  tc^ward  pass  to 
obt  in  the  earliest  occurrence  (Ej)  of  successive  events 
<  using  the  formula 

until  the  end  evenfc(M)  is  reached.  Assuming  that 
the  latest  occurrence  of  the  end  event  (  Lh )  is  the  same 
as  its  earliest  occurrence  (£n),  backward  pass  calcula¬ 
tions  are  made  to  find  the  latest  occurrence  of  each 
predessor  event  until  the  start  event  is  reached  using 
the  following  formula: 

Li 

Thus,  for  each  event,  a  measure  of  slackness  (  SL) 
is  obtained  by  the  difference  (Li-ET^).  Correspondingly, 
the  total  float  (TTij)  in  activities  are  calculated  from 
the  formula  . 

T^J  *  Lj-Ei-dq 

On  this  basis,  the  activities  having  zero  total  float 
are  called  critical  and  a  chain  of  such  activities  from 
the  start  to  the  end  of  the  network  is  called  the  critical 
path.  Several  types  of  float  calculations  may  be  made  for 
various  purposes.  The  critical  path  (or  paths)  provides 
for  management  by  exception  among  other  benefits  derived 
*  from  this  approach. 

Example  (1.1)  There  are  fiva  activities  A,  B,  C,  D 
and  E  having  the  durations  5,  2,  2,  7  and  3  respectively. 


-  y {4 


A  precedes  C  and  D.  B  and  C  precedes  E.  This  Implies 
that  A  and  B  can  commence  simultaneously  while  the  com¬ 
pletion  of  both  C  and  D  terminates  the  project.  The 
following  network  (Fig.  1.1)  Includes  all  calculations 
and  the  critical  path  contains  activities  A  and  D  In 
sequence. 


Fig.  1.1  Critical  path  calculations:  deterministic  case 
In  contrast  to  this,  the  .j.oup  tackling  a  very  large 
project  relating  to  the  Polaris  Missile  had  based  their 
system  on  Event  Orientation.  Since  the  project  was 
highly  complex  and  the  first  of  its  kind,  it  was  not 
possible  to  determine  the  duration  of  individual  activ¬ 
ities  precisely.  The  activity  duration  is  assumed  to 
follow  a  8eta  distribution  .  a-l 

and  that  the  three  time  estimates  t| ,  ta  and  of  (X, 
(optimistic),  the  mode,m0(rnost  likely)  and  b  (pessimistic), 
respectively,  are  available. 


Appealing  to  the  "empherlcal  numerical  manipulation” 
and  truncating  the  normal  distribution  at  +2.66  standard 
deviation,  the  formulaefor  mean  (dl  and  variance  (V) 

areglven  by  r-  .  .-a** 

JL-.  X(t,+-4t2  +  tg)  V  =  [fe  ^  3  J 

The  expected  value  of  the  project  duration  is  cal¬ 
culated  by  Identifying  the  critical  path  as  before  using 
the  expected  activity  durations.  The  distribution  of  the 
project  duration  Is  assumed  to  be  normal  with  the  mean 
value  calculated  as  above.  The  variance  at  the  project 
level  is  a  simple  addition  of  the  activity  variances  ^ 

along  the  critical  path.  (  Mean*  5+7.11  J  " 


Pin.  (1.21:  Critical  path  calculations:  Probabilistic  case 
A  description  of  this  approach  is  given  under  the 
title.  Programme  Evaluation  and  Review  Technique  in  [80], 
the  authors  being  Malcolm,  Roseboom.  Clark  and  Fazer. 
Several  authors  have  examined  the  PERT  assumptions  rather 


critically  and  discussions  are  still  current  as  can  be 
seen  In  Section  Three.  This  type  of  network  is  therefore 
probabilistic  and  event  oriented.  These  two  papers,  l.e., 
Kelly  &  Walker  [64]  and  Malcolm,  et  al  [80],  mark  the 
first  stage  of  development  in  C.P.A.  as  given  in  exhibit  [1], 

1.2  SECOND  STAGE 

The  exhibit  [2]  represents  the  second  stage  of  develop¬ 
ment.  At  this  stage,  importance  is  given  to  the  trade  off 
between  cost  and  time  values  of  individual  activities  in 
such  a  way  as  to  minimize  total  project  costs.  This  sec¬ 
tion  has  been  particularly  attractive  to  operations  re¬ 
search  workers  since,  under  varying  assumptions,  the 
problem  may  be  approached  in  terms  of  one  or  more  well 
known  techniques  of  optimization.  The  two  most  important 
papers  relating  to  this  are  due  to  Kelly  L65]  and  Fulkerson 
[39].  These  two  papers  assume  that  the  cost  vs.  time 
distribution  is  bounded  at  both  ends,  continuous  and  con¬ 
vex.  The  network  flow  approach  of  Fulkerson  also  admits 
all-integer  values  for  which  an  all-integer  solution  is 
obtained.  Several  authors  have  contributed  to  the  ex¬ 
tension  of  this  basic  model.  A  brief  review  of  these 
papers  may  be  seen  in  Section  4  of  this  paper. 

1.3  THIRD  STAGE 

The  second  stage  of  development  providing  a  time- 
cost  trade  off  assumes  an  unlimited  availability  of 
various  resources.  Besides,  the  various  elemental  costs 
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pertaining  to  each  activity  are  not  indicated.  The 
third  stage  of  development  includes  the  time/cost  con¬ 
siderations  in  greater  detail  still  independent  of  re¬ 
source  restrictions.  The  divisibility  of  activities, 
network  decomposition,  use  of  alternative  paths  and 
several  other  refinements  are  shown  in  exhibit  [3]. 

The  time/resource  problem  as  outlined  here  is  basically 
one  of  scheduling  appropriate  resources  to  different 
activities  in  such  a  way  as  to  satisfy  one  or  more  ob¬ 
jectives.  Most  of  the  published  material  on  C.P.A.  are 
included  in  these  three  exhibits,  [1],  [2]  and  [3],  in 
one  form  or  another.  In  particular,  the  exhibit  [3] 
brings  together  recent  publications  in  the  four  areas 
discussed  in  the  next  five  sections  of  this  paper:  Net¬ 
work  algebra,  PERT  assumptions,  cost/time  trade  off  and 
resource  allocation  (heuristic  and  analytical). 

1.4  FOURTH  STAGE 

The  exhibit  [4]  is  an  attempt  to  provide  directions 
of  research  in  C.P.A.  The  i nter-connecti 6ns  shown  in 
the  exhibit  represent  only  a  few  among  several  possible 
model  refinements.  For  example,  the  problem  of  resource 
allocation  has  received  considerable  attention  recently 
(see  sections  5  and  6).  Heuristic  procedures  have  become 
sophisticated  enough  to  include  several  realistic  features 
of  a  project  while  the  available  analytical  models  are 


10 


jStage  I V 1 


EXHIGIT  IV 


% 

i 

£ 


<V  »y„ !  '■I.4  LVtkJ  W  1 


are  unable  to  cope  with  large  complicated  projects. 

However,  attempts  are  being  made  to  combine  costs  and  re¬ 
sources,  using  a  heuristic  as  well  as  an  analytical  frame¬ 
work.  The  source  material  for  the  three  stages  of  develop¬ 
ment  .nd  directions  of  research  as  depicted  in  exhibits 
[ 1 ] »  [2],  [3]  and  [4]  may  be  found  in  these  sections.  For 
example.  "Divisible/Indivisible  activities"  referred  to 
in  exhibit  [3]  is  derived  from  an  article  by  Jewell  [61] 
titled  "Divisible  Activities  in  Critical  Path  Analysis" 
as  reported  in  section  7>  and  indicated  in  exhibit  [5]. 
Similarily,  "Conjunctive/Disjunctive  activities"  of  ex¬ 
hibit  [3]  is  related  to  the  two  papers  by  Roy  and 
Sussmann  [104]  and  Balas  [2]  which  come  under  network 
algebra  (exhibit  [5]) . 


SECTION  II:  NETWORK  ALGEBRA 


2.1  AN  OVERALL  SURVEY 

Irrespective  of  whether  a  network  is  «vent-oriented 
or  activity-oriented,  the  project  network  as  conceived 
by  Kelley  and  Walker  [64]  and  Malcolm,  et  al  [80]  is 
directed,  acyclic  with  activities  having  deterministic 
or  probabilistic  durations.  While  conforming  to  this 
pattern^  Charnes  and  Cooper  [15]  provide  an  interesting 
approach  to  the  critical  path  calculations  using  what 
they  call  a  "directed  subdual  algorithms  One  of  the 
early  papers  suggesting  a  variation  from  the  basic  pat¬ 
tern  is  due  to  Eisner  [31]  who  introduced  decision  boxes 
in  the  letwork.  Elmaghraby  [33]  further  generalized 
this  by  introducing  self  loops  and  feedback  loops. 
Chilcott  and  Thursfield  [17]  illustrate  applications 
of  these  ideas  in  research  management  and  project  evalu¬ 
ation.  Crowston  and  Thompson  [22]  introduce  the  title, 
"Decision  C.P.M."  for  such  networks  containing  decision 
nodes.  Confining  their  attention  to  deterministic  branch 
ing  into  activities,  they  provide  algorithms  employing 
time-cost  trade  off  principles  for  reducing  this  type  to 
a  simple  network  having  unique  activities.  Yet  another 
change  in  the  original  network  is  in  the  use  of  disjunc¬ 
tive  constraints  due  to  Roy  and  Sussmann  [104],  further 
developed  by  Balas  [2],  Even  though  these  disjunctive 
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constraints  are  currently  used  in  simple  machine¬ 
sequencing  problems  illustrated  in  a  network  form, 
the  technique  described  by  Balas  has  much  wider  po¬ 
tential  applications  in  resource  allocation.  This 
section  on  Network  Algebra  (see  exhibit  5)  also  in¬ 
cludes  contributions  by  Parikh  and  Jewell  [93]  and  by 
Luttgen  [78]  on  network  decomposition,  Fisher  and 
Nemhauser  [35]  on  multicycle  project  planning  and 
Schoderbek  and  Digman  [105]  or,  Synthesis  of  PERT  and 
Line  of  Balance.  Even  though  these  papers  relate  to 
the  original  directed  acyctic  network  they  have  a 
significant  bearing  upon  netowrk  topology.  The  follow¬ 
ing  paragraphs  in  this  section  cover  some  of  these- 
articles  in  greater  detail.  Further  discussions  on 
[2,  22,  35  and  93]  are  giver,  in  Section  4. 

2.2  DIRECTED  SUBDUAL  ALGORITHM 

Charnes  and  Cooper  [15]  convert  a  project  graph. 

Fig.  (1.1),  to  a  flow  network  in  which  a  hypothetical 
unit  of  flow  (+1)  leaves  the  start  event  .1  and  enters 
the  end  event  4  .  A  flow  of  (-1)  is  ascribed  to  the 
end  event.  All  intermediate  events  (u^  2  and  3  )  play 

the  role  of  transsii pment  points  having  conservation  of 
flow.  Accordingly,  the  activity  duration,  (  dUj  ),  is 
interpreted  as  the  time  (or  cost)  of  transporting  a  unit 
of  flow  from  event(  i  )  to  event  (j  ).  The  related  linear 
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programming  problem  can  now  be  formulated  with  the  help 
of  non-neg-ti ve  flow  variables  (  Xcj)  corresponding  to 
activities  (  i^j  ),  which  maximizes  the  function 

subject  to  the  terminal  event 
conditions  ]>. Xjj  =  “11  »  I  and  the  flow  conservation 

in  transshipment  points  llAj 3  for  each 

. n-i .  This  generalized  form  assumes  that 

Xij-0  if  there  is  no  activity  connecting  everts  (  t  ) 
and  (j  ) . 

With  reference  to  Fig  («*0  ,  the  L.P.  formulation 
for  the  primal  and  dual  are  as  follows: 

Primal  Problem 

Maximize  f(x)«  5  +  »,  X,j  +  *•  Hi*  ^  *14  +  3  Xi4 

Subject  to 

X,i  +  *13  =  1 

-X|i  +  *13  +  *14  =  ° 

-  x,3  -  X*3  -4*  X34  «  o 

—  XZ4  -  X34  = 

xcj  ^  O,  oJS.  (ij) 

Dual  Problem 

Minimi ze  <jr  (w)  -  -^4 

Subject  to 


UT, 

ur, 


*°2  >  \ 

“^3  >  ^ 

wz  -  ur3  ^  ^ 

10  2  -  UJ4  >  7 

^3  -  >  3 


Where 


—  oO  ^  ^  °0 


~  i;  3  ;  4 
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The  dual  form  has  the  advantage  of  involving  only 
two  dual  variables  in  each  constraint  (teprimal  variable) 
and  therefore  by  assigning  a  value  to  W,  ,  the  absolute 
value  of  can  be  minimized  simply  by  satisfying  the 
dual  constraints  one  by  one. 

Associated  with  each  event  (&),  there  is  a  dual 
variable  (  w* )  whose  value  is  given  by 

Kl-  «} 

where  the  set  (  t, . j  )  represents  the  set  of  all  im¬ 

mediate  predecessor  events  to  event  (k),  Starting  with 
(  iaT|  =  o),  the  values  of  (  U)fc)  ,  can  be  obtained 

by  repeated  use  of  the  above  formula. 

Accordi ngly , 


UJ2 

UJ,  - 

5  , 

wr  = 

-5 

OJ  3 

U7,  - 

w*- 

-7 

IaJ  3 

UJ2  " 

■Z  / 

UJ4. 

" 

-n 

11 

W4 

4 

U3 

-3} 

In  the  usual  critical  path  analysis,  as  given  in 
Section  (1.1) .  (  )  ^  dlcj and  {  L]  -UL  )  ^  dUj 

and  the  critical  path  is  identified  by  a  sequence  of 
activities  (  t,j  )  such  that  (  E»t * L{  ) ,  (  Ej  *  Lj  )  and 
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(£j dij) .  An  analagous  calculation  of  the  dual 
L.P.  problem  would  imply  that  for  those  dual  constraints 

i 

j  which  are  satisfied  as  equalities  it,  Wj,- *y,  -  etc;  , 

i  ^ 

|  e.g.  first,  third  and  fourth  the  corresponding  primal 

variables  Xj2i  X^and  Xj^must  be  non-negative,  where 

I 

l  as  it  must  be  zero  if  the  dual  constraint  is  satisfied 

t 

1  as  an  equality  (43.*  Xj^-- X3/f »  0) .  Assigning  a  value  0  or  1 

y 

to  the  primal  va^able  Xij  and  checking  the  primal  con¬ 
straint  equation  (  X|23 1  ,  *23*0  ,*24*1  )  it  can  be  shown 
that  -UJtf  Xy  given  that  u>,*0  where  |uy£J 

is  interpreted  as  the  minimum  project  duration 
id.  -  (-te)  *»  (5x0+  (2*6)  +  (2  *0)  •+  (-7111)+  (3*0) 

The  logic  given  above  is  contained  in  the  directed 
subdual  algorithm  of  Charnes  and  Cooper,  who  defines,  for 
a  network,  a  spanning  tree  which  forms  a  connected  graph 
with  no  loops.  Accordingly  every  event,  will  have  one 
and  only  one  activity  of  the  project  negatively  incident 
on  it.  Furthermore,  the  spanning  tree  must  have  exactly 
one  chain  leading  from  the  unique  start  event  to  the 
unique  and  event  as  in  Fig. (2-2)  and  a  flow  of  one  unit 
passes  through  the  links  of  this  chain,  while  all  flow 
variables  corresponding  to  other  activities  have  zero 
value.  Then  the  chain  of  links  bearaing  non  zero  flow 
variables  is  the  critical  path  ( <•«•  t  X2.4  ). 

j  Another  illustration  of  this  approach  may  be  seen  in 

|  Ref.  [87]  ,  (pp.  136-139). 


Elmaghraby's  important  paper  [33]  accounts  for  a 
significant  departure  from  the  well  known  directed  acyclic 
networks  having  either  deterministic  or  probabilistic 
activity  durations.  Based  on  the  ideas  of  Eisner  [31] 
regarding  the  probabilistic  branching  of  activities  from 
a  given  event  (as  is  usually  found  in  research  and  de¬ 
velopment  where  the  outcome  of  an  experiment  determines 
further  course  of  action),  Elmaghraby  evolves  an  algebra 
for  generalized  activity  networks. Using  the  fundamental 
probability  laws,  he  provides  formulae  for  reducing  ac¬ 
tivities  in  series,  in  parallal,  in  series -paral 1 al , 
seif  loops  and  feed  back  loops  into  single  activities. 

The  work  of  Elmaghraby  [33]  received  further  impetus 
as  a  result  of  three  recent  papers  [97],  [98]  and  [32], 
the  last  one  being  an  added  contribution  by  the  same 
author.  On  the  basis  of  an  unpublished  doctoral  dis¬ 
sertation  by  Whitehouse  [112]  and  a  RAND  memo  by  Pritsker 
[117],  a  new  term  was  coined  for  the  exploration  of 
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stochastic  activity  networks:  "G.E.R.T."  to  represent  a 
Graphical  Evaluation  and  Review  Technique1. 

With  the  help  of  Signal  Flow  Graphs  [52,  57,  58,  76], 
Pritsker  and  Happ  [97]  have  been  able  to  outlive  an  alge¬ 
bra  for  the  solution  of  networks  having  both  probabilistic 
nodes  and  probabilistic  activity  duration.  Elmaghraby's 
second  article  [32,]  also  includes  these  extensions  in  net¬ 
work  algebra.  The  formalization  of  node  symbols.  Fig. (2. 31) 


Soda  Characteristics  and  Symbol* 


Exclushe-or — The  realization  of  any  branch  leading  into  the  node 
causes  the  node  to  be  realized;  however,  one  and  only 
one  of  the  branches  leading  into  this  node  can  be  real¬ 
ized  at  a  given  time. 

Inclusive -o' — The  realization  of  any  branch  leading  into  the  node 
causes  the  node  to  be  realized.  The  time  of  realization 
the  smallest  of  the  completion  times  of  the  activities 
leading  into  the  Indusne-or  node, 
and— The  node  will  be  realized  only  if  all  the  branches  lead¬ 
ing  into  the  node  are  realized.  The  time  of  realization  is 
the  largest  of  the  completion  times  of  the  activities 
leading  into  the  and  node. 

Deterministic— All  branches  emanating  from  the  node  are  taken  if  the 
node  is  realized,  that  is,  all  branches  emanating  from 
tins  node  have  a  p-parameter  equal  to  one. 

Probabilistic—  At  most  one  branch  emanating  from  the  node  is  taken 
if  the  node  is  realized. 


Fig.  2.31 

and  th2  fundamentals  of  network  reduction  for  deterministic 
activity  durations.  Fig.  (2.32)  are  given  in  these  articles. 

With  regard  to  stochastic  activity  duration,  Pritsker 
and  Happ  [97]  consider  the  principles  of  network  reduction 
for  exclusive  or  nodes  only:  series,  parallal,  series- 
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Original  Net  work 


b 


'  G>-4i>L-<D 


Name  Equivalent  Are 

0-^-0 


Series 

pt  "  pt'Pu 

Parallel 

p.~p(ar\b) 

“And” 

I,  "  max  (/.,  k) 

p,-P{a\Jb) 

Parallel 

/,  ”  [-'.A.  +  hpb  +  1  min  (/,;  fe) 

“Inclusive-or” 

-t.-h\p(ar\b)}/p. 

Parallel 

p.  -  P(aKJb)  -  PiaHb) 

"F.xclusiveor” 

T.  “  [/«/v  +  kPb  + 

M-P(af\b)\/p. 

p.  “  piK  1  ~  p.) 

Self  Loop 
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Fig.  2.32 


paralial,  and  feedback  loops.  Essentially  they  employ 
a  moment  generating  function  (M.6.F.);  Mk(s),  continuous 
or  discrete^associated  with  activity  duration  (t) 
which  is  further  transformed  by  the  conditional  prob¬ 
ability  [p)  oriented  to  the  branch  to  form  a  W-function 
w  (s)  =  f>  Mt  (b)  .  The  outcome  of  this  two  stage  transforma¬ 
tion  results  in  a  system  of  linear  independent  equations 
which  is  amenable  to  fl owgraph-techniques .  The  two  key 
features  of  network  reduction  at  the  transform  level  are 
(i)  The  M-G-F  of  the  sum  of  time  values  is  the  product 
of  tne  M.G.F.  of  individual  time  values  (ii)  the  M-6-F 
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of  a  mixture  of  two  distributions  is  the  sum  of  M-G-F 
of  individual  distributions  each  one  being  weighted 
by  their  conditional  probabilities.  At  the  two  stage 
transformation  level  this  is  reduced  to  simple  addition 
of  the  corresponding  functions  as  in  Fig.  (2.33).  Com¬ 
plicated  feedback  networks  are  reduced  with  the  help 
of  Mason's  Rule  [118,  119]. 
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Network  Type 
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Fig.  2.33 

Pritsker  and  Happ's  article  [97]  also  provides  a 
systematic  procedure  for  converting  an  ' incl usi ve-or '  and 
'and'  type  of  nodes  into  equivalent  ' excl usi ve-or '  types 
given  that  the  activity  duration  is  deterministic.  For 
example,  an  'And'  node.  Fig.  (2.34),  can  be  replaced  with 
two  ' Excl usi ve-or '  nodes,  Fig.  (2.35),  one  representing 
the  eventuality  when  the  'And'  node  is  realized,  the 


Fig.  2.34 


Equivalent  Network  for  on  And  Node. 


Fig.  2.35 

Similarly,  an  ' Inclusi ve-or'  node.  Fig.  2.36,  may  be  re¬ 
placed  by  a  set  of  1 Excl usi ve-or '  nodes.  The  original 
RAND  memo  of  Pritsker  [117]  outlines  possible  approaches 
for  problems  of  this  type  when  the  activity  duration  is 
probabilistic.  Whitehouse  [112]  has  investigated  the 
use  of  tables,  inversion  integrals,  Pearson  curves  and 
other  approaches  for  obtaining  the  distribution  function 
from  an  M-G-F  or  other  exponential  transforms  so  as  to 
ascribe  confidence  limits  to  system  performance.  The 
authors  of  [97]  conclude  that  further  investigation  is 
underway  with  regard  to  the  sensitivity  of  the  system 
performance  measures  to  changes  in  parameters  of  indi¬ 
vidual  branches. 

Fig.  2.36 

A  Simple  Network  with  on  Indcdv-oe  Node. 
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Fig.  2.37 

Equivalent  Network  for  on  IndvWve-or  Node. 

In  a  companion  article  [98],  Pritsker  and  Whitehouse 
present  examples  of  multiple  feedback  loops,  Multiple 
input  and  output  nodes,  selected  problems  in  the  prob¬ 
ability  theory,  a  model  of  a  manufacturing  process  and 
finally  a  problem  in  research  and  development,  originally 
conceived  by  Graham  [46].  In  terms  of  actual  problems, 

GERT  has  been  applied  to  che  analysis  of  a  space  vehicle 
count  down  and  the  zone  refining  of  semi-conductor  material. 

Hornbaker  [56]  applies  GERT  to  the  bond  refunding 
decision  problem  originally  defined  by  Bierman  [5]  as  a 
Markov  Process  where  the  state  of  the  process  is  described 
by  an  interest  rate  in  each  time  period.  Assuming  that 
the  firm  breaks  even  refunding  at  the  current  bond  in¬ 
terest  rate,  Bierman  introduces  a  pay-off  or  reward 
matrix  whose  elements  represent  savings  from  refunding 
each  interest  rate  below  the  break-even  rate  of  interest 
for  each  time  period.  By  converting  the  qualitative 
description  of  the  bond  refunding  process  to  a  model  in 
a  network  form,  Hornbaker  applies  GERT  demonstrating  the 
equivalence  of  results  for  the  same  data  used  in  Bierman's 
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analysis.  Even  though,  for  this  problem,  Bierman's 
simple  approach  is  a  better  one,  GERTholds  promise  in 
the  solution  of  such  problems  should  they  be  made  more 
complex . 

2.4  JOBSHOP  OPERATIONS 

An  entirely  different  point  of  view  of  network 
algebra  is  evident  from  the  works  of  Carol  1  [ «3]  and 
Trilling  [107].  There  are  several  job  shop  operations 
such  as  the  building  of  large  steam  turbines  which  re¬ 
quires  procedures  of  great  complexity.  Trilling  [107] 
describes  a  coding  procedure  based  on  binany  numbers 
while  defining  the  networks  represented  by  the  routings 
or  line  ups  of  shop  orders.  The  combinational  routing 
shows  the  sequence  in  which  various  operations  are  per¬ 
formed  on  simple  parts,  combination  of  parts,  assemblies 
and  disassemblies.  Whenever  fabrication  involves  com¬ 
binational  routings,  one  basic  consideration  is  how  well 
coordinated  the  parts  are  that  come  together  for  assemblies. 
Delay  in  a  part  or  a  sub-group  may  cause  a  delay  in  the 
job  order  as  in  the  critical  path  and  in  some  cases  ag¬ 
gravate  the  work -i n-process  inventory  situation.  Trilling 
consi  *rs  several  decision  rules  in  job  shop  simulations. 

In  particular,  the  works  of  Rowe  [120]  and  Caroll  [13] 
are  significant. 
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2.5  TABLEAU  METHOD 

A  recent  paper  by  Riggs  and  Inoue  [100] 
provides  an  interesting  diversion  from  the  main  theme 
of  network  analysis.  They  suggest  the  use  of  a  tableau 
having  one  row  and  one  column  for  each  activity.  With 
this  format,  the  network  logic  can  be  represented  by  fill¬ 
ing  the  appropriate  cells  in  this  square  tableau,  leaving 
empty  all  the  other  cells  as  in  Fig.  (2.51).  The  various 
time  calculations  are  then  performed  by  extending  this 
table,  introducing  more  columns. 

As  is  evident  from  Fig.  (2.51),  the  tableau  method 
eliminates  the  need  for  representing  separately  the  events. 
The  equivalence  between  event  occurrence  times  and  ac¬ 
tivity  start  or  finish  times  is  established  in  the  note 
of  Fig.  (2.51)  as  well  as  in  Fig.  (2.52)  which  illustrates 
in  the  time  scale,  the  values  of  total  float  and  inde¬ 
pendent  float  for  a  given  activity. 
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cost-time  trade  off  as  well  as  resource  allocation 
problems.  They  depend  heavily  upon  the  time  scale  dia¬ 
gram  rather  than  the  tableau  to  illustrate  the  compu¬ 
tational  features.  However,  the  applicability  of  the 
tableau  method  in  other  areas  is  demonstrated  by  the 
classical  "Stage  Coach  Problems"  which  is  usually  solved 
by  dynamic  programming. 

2.6  DIVISIBILITY  OF  ACTIVITIES 

Jewell  [61]  considers  a  project  where  an  activity 
having  a  total  time  (U)  may  be  devided  and  allocated 
with  non-negative  values  (  t'tj  )  to  two  or  more  locations 
(  C^j  )  denoted  by  the  subset  D,  of  the  project 

network  and  for  each  "t  ,  . N),  represents 

its  time  of  occurrence.  The  original  primal,  and  its 
dual,  as  given  in  Fig.  (2.61)  are  modified  by  bringing 
the  primal  constraint  for  (U)  into  the  objective  function 
with  the  help  of  an  unknown  LaGrange  multiplier  (Q).  The 
modified  form,  Fig.  (2.62),  facilitates  the  use  of  a  stan¬ 
dard  network  flow  approach.  A  parametric  procedure  is 
initiated  with  a  large  value  cf  Q.  As  this  value  is  de¬ 
creased  gradually,  there  is  a  resultant  increase  in  the 
allocation  of  time  (  tcj  )  to  the  D  activities.  The  pro¬ 
cedure  terminates  when  this  allocation  reaches  the  desired 


In  the  case  of  several  divisible  activities  (K), 


the  modified  primal  objective  function  becomes 

The  dual  of  this  may  be  recognized  as  a  multi  commodity 
network  flow  problem.  Jewell  suggests  further  extensions 
of  the  model  to  include  time-cost  trade  off  principles. 
Specific  procedures  to  achieve  these  extensions  are 


still  awaited. 
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SECTION  III:  PERT  ASSUMPTIONS 
3.1  THE  CLASSICAL  CASE 

This  section  is  mainly  devoted  to  directed  acyclic 
networks  whose  activity  durations  are  described  by 
probabilistic  distributions.  The  classical  paper  of 
Malcolm  et  al  [80]  assumes  a  Beta  distribution  for  ac¬ 
tivity  duration.  The  expected  value  and  the  variance 
of  activity  durations  are  derived  from  three  easily 
observable  points  in  the  curve,  the  two  ranges  and  the 
model  value  by  means  of  two  simple  formulas  as  given 
in  Section  (1).  The  econo-iic  considerations  of  informa¬ 
tion  gatherings  and  other  attendant  problems  related  to 
the  complex,  new  and  large  project  in  their  hand  hau 
led  them  to  making  several  approximations  both  at  the 
activity  level  and  at  the  project  level,  i_e  the  manner 
of  aggragating  individual  activity  time  values  through 
the  network  which  lead  to  a  probability  distribution  of 
the  project  duration.  The  authors  then  suggest  that 
the  probability  of  completing  a  project  by  a  given  date 
can  be  computed  by  calculating  the  critical  path  using 
and  expected  activity  duration  as  deterministic  quanti¬ 
ties  and  then  invoking  the  central  limit  theorem. 

This  approach  is  justified  by  the  assumption  that 
a  critical  path  is  comparatively  longer  than  any  other 
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path  in  the  network  so  that  the  realization  of  a  dif¬ 
ferent  path  is  least  probable.  Besides,  the  application 
of  the  central  limit  theorem  is  based  on  the  observation 
that,  in  practice,  the  critical  path  would  have  enough 
number  of  activities  to  justify  the  same. 

3.2  THE  AGGRAVATION  PROBLEM 

While  discussing  the  PERT  assumptions,  Murray  l 90 ] 
states : 

"This  assumption  states  that  the  amount  of  time 
required  to  complete  a  network  or  to  arrive  at 
a  given  event  within  the  network,  has  the  same 
statistical  distribution  as  the  time  required 
to  complete  the  critical  path  leading  into  this 
event.  Or,  to  state  the  assumptions  more  gener¬ 
ally,  the  distribution  of  the  maximum  of  a  set 
of  normally  distributed  random  variables  is  equal 
to  the  distribution  of  the  random  variables  with 
the  larges  first  moment. 

This  assumption  was  incorporated  into  PERT 
because  it  was  needed  to  make  possible  the  cal¬ 
culation  of  the  expected  completion  times  for 
varioas  events  in  the  PERT  networks.  This  as¬ 
sumption  was  known  to  be  invalid  and  to  introduce 
a  bias  in  the  event  completion  date  calculations 
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but  it  was  thought  that  the  biases  were 
small.  Work  has  been  done  towards  ulti¬ 
mately  correcting  this  assumption.  Clark 
[1*1]  has  developed  approximate  methods 
for  computing  the  moments  of  the  maximum 
of  a  set  of  normally  distributed  random 
vari abl es . " 

Given  a  simple  network  consisting  of  only  two 
parallal  activities  having  durations  (  t|  )  and  (  } , 

means  (  pi  )  and  (  pa,)  and  variances  (  6r?x)  and  {  ) , 

Clark  [19]  shows  that  for  the  random  variable  that  is 
the  maximum  of  (  tt  )  and  (  ,  the  mean  (  pM)  and 

the  variance  )  are  given  by 

|*M  *  Hi  4>(*0  4  4  C"*)4  O-  & 

and  «  >.  ii^ 
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For  the  case  of  equal  means  ,  as  in  Tippet  [106], 

JJL,  »  p- 2  ,  -  ~  0  <Vyv*“ 

(A-m  =  +  (:r^)J6",'+625'  ‘kr'J- 

«rf-  - 

so  even  for  this  simple  case,  the  mean  and  variance  of 
the  distribution  of  the  maximum  are  dependent  upon  the 
variances  of  the  two  activities. 

Clark  [19]  then  assumes  that  the  elements  in  the 
network  are  normal  random  variables,  that  the  paths 
through  the  network  are  correlated,  that  the  paths  have 
unequal  means  and  variances  and  that  the  distribution 
of  the  larges*  value  is  approximately  normal.  For 
such  a  case,  he  gives  an  iterative  procedure  whereby 
one  integrates  his  way  through  the  network  eliminating 
an  element  at  a  time  in  terms  of  equivalent  elements 
and  finally  ends  up  with  the  expected  value  and  variances 
of  the  network.  While  summarizing  Clark's  bias  correction 
procedure  as  given  above,  Klingel  [70]  observes  that  his 
approach  is  very  tedious  and  also  an  approximation  since 
it  is  necessary  to  assume  normality  at  each  iteration. 
Besides,  the  error  of  the  approximation  has  not  yet  been 
clearly  established.  Moder  and  Phillips  [87]  (pp..  229- 
239)  provide  an  illustrative  application  of  this  procedure. 
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3.3  COEFFICIENT  OF  SKEW 

In  support  of  the  original  paper  by  Malcolm  et  al 
[80],  Clark,  one  of  the  authors  of  this  paper,  provides 
a  practical  explanation  in  [20],  besides  an  attempt  to  validate 
the  probability  statements  in  terms  of  a  truncated  normal 
distribution  as  an  appropriate  simple  model  for  specifying 
che  ratio  of  the  standard  deviation  to  the  range.  However, 
Grubbs  [47]  has  pointed  out  the  subjective  nature  of  the  PERT 
estimation  problem  and  the  restrictions  on  the  Beta  Distri¬ 
butions  that  are  required  to  validate  the  use  of  PERT  means 
and  variances. 

Referring  to  the  PERT  formulae  given  in  the  first 
secti on-Grubbs  shows  that  She  parameters  (f>)  and  ty) 
of  Beta  distribution  are  restricted  to  the  values. 

f>--  3  ta  ,  <£=  3  +  fz  o<  + 

In  particular,  he  indicates  the  existence  of  a  coefficient 
of  skew  M  given  by  the  following  equation. 

y,  »  Ms/^3  =■ 

This  is  indeed  inconsistent  with  the  statement  by 
the  original  authors  that  no  assumption  is  made  about 
the  position  of  mode  relative  to  the  extreme  values. 

Donaldson  [29]  suggests  the  use  of  the  mean  (  ft\0  a  tj,  ) 
rather  than  the  mode  (  m  )  in  addition  to  the  two  limits 
of  the  distribution  {  a  »  fc,  and  b=  £3  )  as  the  initial 
estimates.  Accordingly  rt)0^  4-  /(£+  and  V7 

(variance  estimator)  =  D/o>+^o 

Assuming  that  Jp>2/  and  <£> Z  ,  the  smallest  value  of 
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(  p  4  )  which  satisfies  the  necessary  condition*  ic 

2  +  3K+  S')  »  where  S  is  a  small  positive  value  and 
K  *  r***  (t^-tO/Ctj  -  'Q 

As  pointed  out  by  Coon,  a  referee  of  this  paper,  t..e 
family  of  curves  generated  by  this  approach  do  not  lead  to 
complete  generalization  of  PERT  activity  time  distributions, 
even  though  there  is  a  considerable  relaxation  of  the  current 
restrictions. 
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3.4  BIASES  DUE  TO  ACTIVITY  SUBDIVISION 


In  another  early  paper  on  this  topic,  Healy  [54] 
investigates  the  bias  introduced  in  the  probability 
statements  regarding  project  completion  if  the  network 
activities  are  subdivided  as  from  Fig.  ( 3  •  A-l )  "to 
&g(3'4^).  The  completed  probabilities  vary  depending 
upon  whether  the  scheduled  date  for  an  event  is  later 
or  earlier  in  time  than  the  expected  date.  An  important 
reason  for  this,  as  Healy  points  out,  is  that  the  stan¬ 
dard  deviation  of  the  sum  of  two  random  variables  (square 
root  of  the  sum  of  two  variances)  is  not  the  same  as  the 
standard  deviations  of  the  combined  distribution.  If 
the  basis  of  the  calculation  as  in  PERT  is  based  upon 
the  approximation  that  for  unimodal  distributions,  stan¬ 
dard  deviations  can  be  estimated  roughly  as  one  ninth  of 
the  range. 

Even  though  this  represents  an  interesting  mathe¬ 
matical  excursion  into  probability  calculations,  Clark 
[54],  a  referee  cf  this  paper,  comments  that  the  decision 
to  subdivide  an  activity  is  determined  by  the  nature  of 
the  activity  and  thac  the  decision  is  not  a  matter  of 
analysis.  Admittedly,  it  is  not  always  easy  to  deter¬ 
mine  the  appropriate  amount  of  subdivision.  Clark  argues 
that  this  determination  should  be  referred  to  the  infor¬ 
mation  source,  not  the  analyst.  Millstein  [54],  another 
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referee,  comments  that  these  probability  estimates 
should  be  compiled  with  other  factors  such  as  the  vari¬ 
ances  for  expected  event  time,  latest  event  time,  schedule 
event  time,  the  amount  of  event  slack  and  the  density  of 
activities  leading  to  the  event.  He  further  observes 
that  in  normal  operating  applications  of  PERT,  an  in¬ 
crease  in  the  density  of  activities  from  3  to  6  times 
the  original  planned  schedule  is  abnormal. 

3.5  i'OINT  DISTRIBUTION  OF  IMMEDIATE  PREDECESSOR  ACTIVITIES 
Another  important  early  contribution  to  PERT  as¬ 
sumptions  is  due  to  Fulkerson  [40].  One  of  his  signi¬ 
ficant  observations  in  this  direction  is  that  the  PERT 
estimate  (<J)  (ie^  use  of  expected  activity  duration  in 
the  deterministic  calculation  of  the  critical  path  for 
project  duration)  grossly  underestimates  the  true  expected 
project  duration  ( <2. )  which  should  include  the  whole  ac¬ 
tivity  distribution  rather  than  its  expected  value.  Ob¬ 
viously,  it  would  be  impossible  to  calculate  the  true 
mean  for  large  networks  having  continuous  activity  dis¬ 
tributions.  However,  Fulkerson  suggests  a  method  of 
obtaining  a  fairly  good  lower  approximation  [f)  for 
expected  project  duration  for  discrete  activity  dis¬ 
tribution  such  that<j^f^6.  The  basic  assumptions  in 
his  model  are  that  the  jobs  which  immediately  precede 
any  given  job  may  have  a  joint  distribution  of  times 


38. 


p(LB^)  while  there  is  independence  between  jobs  which 

* 

immediately  precede  different  jobs.  For  each  event  t-  , 

(given  t<j)  there  is  a  bundle  of  all  preceding  activities 

(  (5{,) ,  such  that  a  set  of  activity  durations 

"•  can  be  chosen  with  the  appropriate  joint  probability 

distribution  jp(tsi)  en  activity  duration  within  a  bundle* 

Even  though  all  preceding  activities  are  included  for 

generality,  only  the  immediate  predecessor  activities  and 

the  respective  predecessor  events  to  event  (  )  are  used 

in  the  calculations  of  the  expected  occurence  (  jfi  )  of 

#  • 
event  (L)  (^e  length  of  the  path  from  origin  to  event  (C)). 

Fulkerson  defines  these  numbers  (  )  recansively  as 

follows:  f,  =  0 


He,  thens  illustrates  the  computational  features 
through  simple  examples  as  given  below,  proceeding  from 
one  node  to  the  next  successor  node  taking  into  account 
all  discrete  probability  values  for  each  activity  under 
consideration.  When  once  the  expected  length  of  one 
node  is  calculated  this  is  used  as  a  basis  for  the  next 
successor  node.  This  process  is  comparable  to  the  de¬ 
terministic  calculation  of  the  critical  path.  This 
approach  often  results  in  an  improved  estimate  even 
though  more  calculations  are  involved. 
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Example  3.51 :  Let  the  network  be  that  of  Fig.  (3.51) 
and  suppose  that  each  arc  independently  assumes  the 
length  0  or  1  with  probability  1/2  (Ref. [40],  p.  812) 


Fig.  3.51 

The  values  of  (  ) »  *•  3, 4  are  obtained  by  using 

the  expected  activity  duration  (  Acj  =  Vi.  )  all  (t,  j}  )  i_e 

®!,  =  0  )  9 1  -  9,  +  4“'k  y  83  =  ***  *9z4*)  - 1 

and  m*.*  (qa+-L  ,  q3-rr)  =.  X 

The  values  of  ( ^  ),  {,*1,2.,$,+  are  obtained  by 
using  the  recursive  formula  given  above. 
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Thus  the  statement  "9;^  l'z>3  *nJ’4  "  Is 

*  L  ^ 

validated  by  this  example  since  <3=  (.  O ,  ‘/a. ,  l  ,  V*") 
while  f  s  (  O,  ‘/a. ;  % ,  55/3i)  •  A  formal  proof  for 

may  be  seen  in  [403  - 

Clingen  [21J  describes  an  extension  of  Fulkerson's 
algorithm  to  include  the  continuous  case  involving 
Stieltjes  integration.  The  original  formulation  of 
the  continuous  case  is  not  computationally  feasible 
since  in  general  it  would  involvenulti-dimensional  nu¬ 
merical  integration.  However,  he  derives  an  equivalent 
expression  which  for  a  discrete  case  is  reduced  to 
simple  finite  sums  and  for  the  continuous  case,  a  one 
dimensional  integration. 

In  continuation  of  the  last  two  papers,  Elrnaghraby 
[116]  suggests  two  methods  both  of  which  further  improve 
upon  Fulkerson's  estimates  of  the  true  mean  project  dura¬ 
tion.  His  first  method  is  based  upon  the  observation 
that  even  though  the  true  mean  of  project  duration  is 
not  affected  whether  the  calculations  are  made  in  for¬ 
ward  or  backward  directions,  the  intermediate  events 
would  have  different  values  using  Fulkerson's  algorithm. 
Therefore,  by  obtaining  the  two  estimates  of  expected 
length  of  a  subnetwork  from  forward  and  backward  calcu¬ 
lations,  the  higher  of  the  two  values  should  be  chosen 
for  subsequent  calculations.  This  would  eventually  lead 
to  a  better  estimate  of  the  project  duration.  The 
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second  approach  recommenced  by  Elmghraby  is  a  generali¬ 
zation  of  Fulkerson's  procedure,  by  increasing  the  num¬ 
ber  of  point  estimates  for  each  event  and  superimposing 
his  first  method  of  forward  and  backward  probability 
calculations . 

Elmaghraby's  second  approach  would  give  a  still 
better  estimate  of  the  true  mean  project  duration. 
However,  each  improvement  of  original  PERT  calculations 
involves  extra  computation.  Elmagrnaby:s  article  pro¬ 
vides  a  basis  for  an  effective  trade  off  between  compu¬ 
tational  cost  and  accuracy  in  estimating  project  dura¬ 
tion,  since  his  second  approach  can  be  generalized 
theoretically  to  any  number  of  point  estimates  per¬ 
taining  to  each  event,  until  the  probability  distribu¬ 
tion  function  as  a  whole  is  taken  into  account. 

3. 6  INTERMEDIATE  QUANTILES  ftND  GAMMA  DISTRIBUTION 

Fulkerson  [40],  Clingen  [21]  and  tlmaghrapby  [116] 
as  can  be  seen  in  Exhibit  [6]  form  a  stream  of  ideas 
brought  up  to  date.  On  the  other  hand,  Murray  [90]  is 
mainly  concerned  with  the  use  of  the  Beta  distribution, 
three  time  estimates  and  the  approximate  formula  for  the 
standard  deviation.  He  shows  the  dependency  of  the 
formula  for  standard  deviation  upon  the  definition  of 
the  closed  range,  the  optimistic  and  pessimistic  values, 
and  the  unimodal  beta  distribution.  Rather  than  having 


defined  boundaries,  he  suggests  estimation  of  inter¬ 
mediate  quantiles.  He  claims  that  this  revision  would 
more  fully  utilize  the  experience  of  the  estimator  and 
at  the  same  time  the  derived  estimates  are  more  easily 
related  to  the  moments  of  the  distribution  under  con¬ 
sideration.  Alternatively,  he  recommends  the  use  of 
the  Gamma  distribution  instead  of  the  Beta  distribution, 
since  the  former  has  the  advantage  of  requiring  only  two 
time  estimates  to  specify  its  parameter.  However,  the 
major  defect  of  the  Gamma  distribution  here  is  that  its 
mode  is  always  less  than  its  mean.  He  has  also  advocated 
the  use  of  probability  papers  in  obtaining  a  rapid  esti¬ 
mate  of  the  bias  introduced  by  PERT. 

3.7  ERROR  ANALYSIS  AT  ACTIVITY  AND  PROJECT  LEVEL 

MacCrimmon  and  Ryarac's  paper  [79]  represents  an  im¬ 
portant  contribution  in  the  systematic  analysis  of  PERT 
assumptions.  The  bias  problem  is  divided  broadly  into 
two  categories:  The  activity  level,  the  project  level. 

At  the  level  of  individual  activities  they  consider  three 
possible  sources  of  error.  (1)  The  Beta  distribution 
(2)  the  standard  deviation  and  mean  (3)  imprecise  three 
time  estimates.  They  derive  analytical  expressions  in 
each  case  to  indicate  the  magnitude  and  direction  of 
error.  For  example,  the  Beta  distribution  is  compared 
with  the  quasi-uniform  distribution  at  one  extreme  and 


the  quasi-delta  function  at  the  other  extreme  to  obtain 
the  bounds  of  error  introduced  in  calculating  the  activity 
mean  (error  =(l/3)(l-2m)  and  standard  deviation  ( */6) . 
Commenting  on  this  section  of  the  paper,  Lukaszewicz  [77] 
suggests  the  use  of  a  quasi-uniform  distribution  of  a 
slightly  different  interval  which  improves  upon  the 
bounds  of  error.  The  second  source  of  error  due  to  sim¬ 
plified  formulae  for  the  standard  deviation  and  mean  are 
obtained  by  comparing  these  estimates  with  actual  values 
of  the  standar  deviation  and  mean  expressed  as  a  func¬ 
tion  of  the  mode  and  oC  ,  one  of  the  parameters  of  the 
Beta  distribution.  The  worst  absolute  error  in  the  mean 

is  J-  (  4»n+l)  —  W (<*+»)/  (o<+2  Vn)  I 
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The  worst  absolute  error  in  the  standard  deviation  is 


—  -  J  yY\2-(c(+  f)  -*-t n)  j  +  2rY5)a'(«rf-»-3Tfc) 

Similarly  for  the  third  case,  the  lack  of  precision  in 

the  estimating  procedure  is  analyzed  by  allowing  the 

three  time  estimates  to  vary  within  reasonable  bounds, 

and  the  effect  of  this  in  the  calculation  of  activity 

mean  and  standard  deviation  is  studied.  The  worst  ab- 
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ile  that  of  the  standard  deviation  i:  Q/30) * 
The  second  part  of  MacCrimmon  and  Ryavac's  paper  is 


solute  error  in  the  mea 
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concerned  with  the  calculations  underlying  the  proje:t 


mean  and  standard  deviation  as  well  as  the  probability 
statements  which  follow  these  calculations.  As  observed 
by  Fulkerson  [40]  earlier,  the  project  duration  is  always 
biased  optimistically  by  PERT  calculations,  while  the 
standard  deviation  may  be  biased  in  either  direction. 

They  indicate  that  the  network  configuration  Influences 
the  magnitude  of  errors  at  this  level.  In  particular, 
the  length  of  the  critical  path  as  compered  to  other 
paths,  the  number  of  sub-critical  paths  very  close  to 
the  critical  one,  the  number  of  common  nodes  among  vari¬ 
ous  patns,  etc.  tend  to  either  compensate  or  aggravate 
the  error  values.  However,  it  should  be  possible  to 
eliminate  certain  activities  or  parallal  paths  from  the 
network  if  the  sum  of  the  maximum  time  along  this  path 
is  less  than  the  sum  of  the  minimum  times  along  another 
existing  path.  In  fact,  this  principle  has  been  ex¬ 
ploited  by  Handa  [51]  while  considering  time/cost  trade 
off  problems,  so  as  to  reduce  the  total*  computational 
effort. 

3.8  USE  OF  MONTE  CARLO  METHOD 

While  MacCrimmon  and  Ryavac  o?  Rand  Corporation 
were  working  cna  comprehensive  analysis  of  PERT  assump¬ 
tions,  Va.i  Slyke  [108],  also  of  Rand  Corporation,  was 
exploring  the  use  of  Monte  Carlo  methods  to  yield  solu¬ 
tions  to  the  Pert  problem.  This  method  has  a  consider¬ 
able  flexibility  in  the  use  of  any  theoretical  or  empirical. 


continuous  or  discrete  distribution  for  activity  dura¬ 
tion.  The  effect  cf  introducing  or  neglecting  highly 
arbitrary  assumptions  relating  to  these  distributions 
can  be  studied.  As  in  other  problems,  the  Monte  Carlo 
approach  here  involves  obtaining  a  sample  value  fo*“ 
each  activity  duration  in  tno  network  and  proceeding 
with  the  usual  deterministic  critical  path  calcul  .ons 
to  oDtain  the  project  duration.  A  large  set  of  trials 
and  further  statistical  computations  would  yield  dis¬ 
tributions  and  parameters  of  interest.  Either  as  a 
means  of  validating  certain  assumptions  or  as  a  basic 
method  of  computing  desired  results,  Monte  Carlo  offers 
unique  opportunities.  Van  Slyke  observes  that  the 
Monte  Carlo  estimate  of  the  mean  project  length  is 
unbiased  and  the  correlation  between  activities  can 
be  efficiently  handled.  Another  important  by-product 
of  this  approach  is  the  availability  of  a  ’criticality 
index'  Tor  each  activity  expressed  in  terms  of  the 
probability  that  the  activity  will  be  on  the  critical 
path.  Van  Slyke  proposes  two  methods  for  computational 
short  cuts,  one  of  which  has  also  appeared  ir  MacCrimmon 
and  Ryavac's  paper,  leading  to  the  elimination  of  cer¬ 
tain  paths,  if  the  activity  duration  has  finite  /anges. 
Fcr  other  types  of  activities.  Van  Slyke  suggests  the 
elimination  of  all  activities  that  were  never  critical 
withi;.  a  limited  number  of  trials. 


3.9  HERE  AND  NOW  FORMULATION 


The  article  by  Jewell  [GO]  provides  an  interesting 
diversion  from  ’.he  type  discussed  in  this  section.  He 
makes  use  of  the  principles  involved  in  the  cost/time 

trade  off  in  the  face  of  uncertain  activity  duration. 

\ 

At  the  outset,  forJ each  activity  a  preliminary  level  of 
effort  measured  in  terms  of  money,  manpower  or  other  re 
sources  is  determined  and  a  probability  distribution 
for  this  level  of  effort  is  given.  With  this  format 
for  each  activity  the  planner  estimates  the  interval 
needed  for  the  completion  of  each  activity.  These  es¬ 
timates  are  then  used  in  the  usual  way  for  the  calcu¬ 
lation  of  the  project  duration.  The  event  times  cal 
culated  at  this  state  are  fixed  once  and  for  all.  In 
Jewell's  "Here  and  Now"  formulation  of  the  problem,  the 
activity  duration  becomes  known  only  when  it  gets  under 
way.  Based  upon  the  difference  between  the  allotted 
time  interval  and  the  actual  duration  of  the  job,  ad¬ 
ditional  effort  may  be  expended  to  stay  on  schedule. 
Assuming  a  convex,  continuous  time-resource  curve  having 
a  quadratic  form  with  defined  boundaries,  he  employs 
the  quadratic  programming  algorithm.  This  leads  to  a 
project  resource  curve  where  the  minimal  value  of  av¬ 
erage  additional  effort  expended  is  plotted  against 
the  project  duration.  The  curves  assume  a  piece  wise 
quadratic  form 


3.10  POLYNOMIAL  APPROXIMATIONS 

An  excellent  theoretical  contribution  in  PERT  cal¬ 
culation  is  due  to  Martin  [82].  His  basic  assumptions 
regarding  the  probability  distribution  of  activity  dura¬ 
tion  are  that  they  are  independent,  have  a  finite  range 
and  are  described  by  a  class  of  polynomial  density  func¬ 
tions.  The  last  assumption,  which  is  the  unique  feature 
of  this  paper  was  chosen  mainly  to  facilitate  digital 
computation.  The  polynomial  approximation  does  produce 
an  error  in  the  calculation  of  project  duration,  but 
Martin's  approach  provides  a  basis  for  a  systematic  study 
of  network  configuration.  In  this  paper,  he  presents  a 
detailed  examination  of  the  convolution  operation  for 
polynomials.  By  performing  all  computations  in  symbolic 
or  algebraic  form,  the  reduction  algorithm  described 
here  circumvents  the  problem  of  trucation  error  usually 
found  in  numerical  computation.  His  reduction  algorithm 
converts  a  series-parallal  network  to  an  equivalent 
single  activity.  Since  his  algorithm  also  includes  in 
a  modified  form  the  transformation  of  any  directed  acy¬ 
clic  network  to  series-parallal  form,  an  efficient  prr - 
cedure  now  availab'e  for  finding  the  distribution  of 
project  completion  time  and  the  probabil i ty  that  an  ac¬ 
tivity  is  on  the  critical  path,  "criticality  index",  as 
originally  conceived  by  Van  Slyke  [108].  The  main  dif¬ 
ficulty  in  the  practical  imp!,  entation  of  this  approach 
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is  in  the  computational  requirements,  in  particular,  the 
exponential  growth  of  polynomial  coefficients  when  the 
number  of  activities  in  the  network  increases  ^53]].  The 
polynomial  coefficients  of  the  completion  times  for 
certain  sub-networks  are  obtained  from  the  recurrence 
formulae  derived  in  his  papers. 

3.11  A  NEW  SYSTEM  OF  NETWORK  CLASSIFICATION 

Yet  another  attempt  to  remove  the  biases  in  PERT  as¬ 
sumptions  is  made  by  Hartley  and  Wortham  (]53j  .  They  pre¬ 
sent  a  new  system  of  network  classification  which  dif¬ 
ferentiates  between  'uncrossed',  'crossed'  and  'multiple 
crossed'  networks.  Their  new  definitions  include  terms 
such  as  an  universal  point,  consecutive  points,  sets  of 
second  ordered  branches  and  so  on.  This  classification 
allows  them  to  adopt  different  solution  procedures  for 
different  types  of  networks  depending  on  their  degree 
of  involvement  and  complexity.  This  method  of  analysis 
includes  an  operational  calculus  such  as  numerical  inte¬ 
gration  for  basic  types  of  networks,  a  Monte  Carlo  pro¬ 
cedure  for  more  involved  networks  and  analytical  solu¬ 
tions  for  simple  networks.  The  novel  feature  of  this 
article  is  the  attempt  to  synthesize  various  contribu¬ 
tions  in  this  direction  and  evolve  a  statistical  theory 
for  the  derivation  of  unbiased  distribution  of  the  pro¬ 
ject  completion  times,  with  a  provision  for  sensitivity 
analysis  relating  to  assumptions. 


... 


3.12  REDUCTION  ALGORITHMS 


In  his  doctoral  dissertation  under  Hartley,  Ringer 
[101]  presents  algorithms  to  reduce  the  original  PERT 
network  to  an  equivalent  network  consisting  of  fewer 
activities.  The  distribution  of  completion  times  for 
K  activities  in  series  are  reduced  to  a  single  distri¬ 
bution  of  the  sum  of  K  random  variables.  If  the  K  ac¬ 
tivities  are  in  parallal,  the  distribution  of  the  maxi¬ 
mum  of  K  random  variables  is  chosen.  Certain  types  of 
crossed  subnetworks  can  also  be  reduced.  These  algorithms 
are  repetedly  applied  until  no  further  reduction  is  pos¬ 
sible  either  because  a  single  activity  remains  or  because 
crossed  networks  a  non  reducible  form  remain.  In  the 
later  case.  Ringer  outlines  an  efficient  stratified 
sampling  scheme  as  a  part  of  the  Monte  Carlo  technique 
to  obtain  the  mean,  variance  and  cumulative  distribution 
function  of  the  completion  time.  In  addition  to  the  above, 
he  also  presents  a  model  which  has  dependency  among  acti¬ 
vities  in  series  or  parallal. 

3.13  EXPERIMENTS  ON  MULTIPLE  PARALLAL  PATHS 

More  recently,  Klingel  [70]  reports  an  interesting 
study  of  a  real  network  whi ch  indicates  the  bias  introduced 
in  project  duration.  The  network  chosen  by  him  has  multi¬ 
ple  parallal  paths  with  relatively  equal  means  and  large 
variances.  Some  of  Klingel's  observations  regarding  the 
works  of  T  ippett  and  Clark  in  the  analytical  estimation 
of  project  duration  have  already  been  locked  into  in  the 
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beginning  of  this  section.  His  article  gives  details  of 
experimental  design  and  simulation  procedures  adopted 
for  a  large  network  of  a  thousand  odd  activities  relat¬ 
ing  to  multiple  restaurant-service  installations  between 
1  and  10,  The  network  configuration  is  changed  in  terms 
of  parallel  paths,  retaining  the  common  constraints. 

With  regard  to  the  nature  of  activity  distribution,  Klingel 
has  chosen  Fisher's  coefficient  of  variation  (the  ratio  of 

a  4 

stji^rd  deviation  to  the  mean)  as  an  appropriate  measure 
of  this  variability.  In  all,  24  combinations  between  the 
number  of  parallel  paths  and  the  Fisher's  coefficient  of 
variation  has  been  looked  into,  making  25  replications 
for  each  combination.  The  final  result  of  this  real  ex¬ 
ample  shows  the  incremental  accumulation  of  bias  in  project 
duration  up  to  50%,  as  the  variability  in  activity  duration 
apH  the  multiplicity  of  paths  are  increased  as  in  Fig. 
(3.13).  In  Conclusion,  Klingel  makes  an  important  obser¬ 
vation  pertaining  to  resource  allocation  in  projects. 

"There  are  a  large  number  of  resource  al- 
*  ^cation  programmes  available  which  reduce  de¬ 
terministic  completion  time  by  increasing  the 
parallelism  of  the  network.  While  increasing 
parallelism  reduces  deterministic  completion 
time,  it  also  reduces  the  probability  of  com¬ 
pleting  on  or  before  the  new  deterministic 
date." 
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Fig.  (3.13) 


ESTIMATING  BEHAVIOR  OF  INDIVIDUALS 

The  two  most  recent  papers  in  this  topic.  King  and 
Wilson  [67]  end  King,  Vittebrongel  and  Hazel  [68]  repre¬ 
sent  another  interesting  diversion  from  the  basic  theme. 
The  second  paper  referred  to  above  is  more  or  less  a 
continuation  of  the  first  one,  validating  the  tentative 
conclusions  reached  earlier.  This  series  fulfills  a 
Song  felt  need  to  initiate  research  on  the  estimating 
behavior  of  individuals  in  rela«ion  to  PERT  assumptions. 
The  research  group  has  collected  a  wealth  of  data  about 
past  subjective  judgements  of  estimators  in  testing  two 
hypotheses  *. 


(1)  Preactivity  time  estimating  accuracy  improves 
as  the  beginning  date  of  the  activity  approaches. 

(2)  Time  estimates  made  during  the  progress  of  an 
activity  improve  in  relative  accuracy  as  the  completion 
date  of  the  activity  approaches. 

Contrary  to  the  general  beliefs  as  stated  in  the  two 
hypotheses,  the  researcher  found  that  the  intuitive  learn¬ 
ing  concept  inherent  in  first  hypothesis  was  totally  lack¬ 
ing  in  the  pre-activity  estimating  behavior  of  individuals 
even  if  they  were  given  ample  opportunities  to  re-escimate 
The  second  paper  reports  the  findings  of  five  different 
government  contracts.  Strangely  enough,  even  the  second 
hypothesis  could  not  be  validated  by  their  test  results. 
Hence  they  conclude  that  there  is  no  significant  change 
in  the  accuracy  of  estimating  the  remaining  portion  of 
an  activity  as  the  portion  of  activity  remaining  becomes 
smal 1 er . 

On  the  basis  of  this  information,  *:he  study  group 
proposes  the  establishment  of  a  super  network  which  in¬ 
corporates  adjusted  time  estimates,  under  the  assumption 
that  there  is  a  basic  stability  in  the  structure  of  the 
estimating  system.  King  and  Wilson  consider  several  cost 
models  in  establishing  the  adjustment  procedure  and  they 
conclude  that  models  based  on  the  minimization  of  total 
error  cost  may  be  useful  in  adjusting  time  estimates  if 
valid  cost  relationships  can  be  established. 


3.15  CHANCE  CONSTRAINED  PROGRAMMING 


The  direct  and  the  dual  linear  programming  approach 
to  the  deterministic  network  time  calculations  are  given 
in  subsection  (2.2)  where  Charnes  ana  Cooper  [15]  pre¬ 
sent  a  directed  sub  dual  algorithm.  The  probler  may 

now  be  rephrased  as  follows: 
n 

Primal:  max  fcj  Xj 

SuflTECT  To 

£.€yVfli  ;  (XJ  >o) 

J*'  (i-o,1, . ra) 

where  &0-  ~l  0  2  I  and  all  other  f  tjis  the  duration 

of  activity  (j  )  and  are  the  incidence  numbers  for 

the  networx. 

i n 

Dual:  min.  (X^ 

t  *  0 

SueiecT  to  7w 

t- Ui  £y  ^ 

C*  o 

Charnes,  Cooper  and  Thompson  [16]  present  a  stochastic 

version  of  this  problem  by  assuming  that  the  tj  values  are 

random  variables.  Accordingly  the  dual  problem  consists  of. 
m 

SUBJECT  to 

l -O 

p  ( £,  u-t  ^  tj)  >  pj 

V  C-*0 

where  P  means  'probability  and  O  lis  a  preassigned 
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« 

probability  measure  for  J  =  \,Z*, --*•  H  .  This  programming 
problem  is  therefore  chance  constrained.  Let  be  the 
marginal  distribution  of  tj  such  that 

P(lui£c j  >/tj)  =  F] 

\  is  ft  '  L=0 


This  leads  to  the  inequality 
"m 


>  FZj  ( £j)  ,0*'^ . *) 

C  =  0 

where  ^  (fc)  represents  the  Fractile  associated  with 
fi j  .  On  this  basis,  the  primal  has  the  form: 


5u67ecx  to 


^£ljXj  -  ftc  ,(Xi»o),fc.o(...)5) 


Using  this  approach,  Charnes,  Cooper  and  Thompson 
[16]  show  the  method  determining  the  critical  path  for 
every  possible  collection*#  t  j  )  realizations  with  their 
associated  probabilities.  Curves,  based  on  electronic 

i 

computation  for  a  simple  network  with  exponential  ac¬ 
tivity  durations  are  also  given.  The  interacting 
properties  of  the  network  show  that  the  overall  be¬ 
havior  of  the  project  duration  is  very  different  from 
that  of  individual  activities,  even  for  this  monotcui cal ly 


decreasing  curve.  The  use  of  this  approach  for  other 
types  of  curves  and  several  other  C-xtenUons  of  the  basic 
model  are  still  under  progress. 


SECTION  IV:  COST/TIME  TRADE  OFF 


4.1  THE  PROBLEM 

Let  the  duration  (  Vcj  )  of  each  activity  (  i,j  ) 
vary  between  the  crash  time  (  (kij  )  and  the  normal  time 
(  bcj  )»  (1_£  4  ^Lj  btj  )  and  1&t  each  unit  de¬ 

crement  in  time  requires  an  additional  cost  (Ccj  >. 
the  cost  slope.  The  activity  cost  curve  is  therefore 
linear,  continuous  and  bounded  at  both  ends.  This  vari 
ability  of  duration  at  the  activity  level  introduces 
several  feasible  schedules  at  the  project  level  for 
each  feasible  project  duration  (  A  )• 

Associated  with  each  schedule  are  the  direct  and 
indirect  project  costs  which  lead  to  total  project  cost 
curve.  Several  approaches  are  available  in  obtaining 
the  direct  project  cost  curve.  This  section  is  related 
to  the  computational  details  of  this  problem  under  vary 
ing  assumptions. 

4.2  USING  THE  CRITICAL  PATH 

A  simple  procedure  involving  the  critical  path 
calculations  is  given  below.  The  example  in  subsection 
(1.1)  is  suitably  modified  to  include  cost  slopes  and 
crash  durations  as  given  in  fig.  (4.10). 


(!)  For  each  activity,  let  3  ^Cj  and  compute  the 
project  duration  (eg/  12  weeks)  and  identify  the 
critical  path  ( ejg_;  1,  2,  4). 

(11)  Choose  among  the  critical  activities,  the  one  with 
minimum  cost  slope  [eg/  min.  (10,  20)  =  10]. 

(ill)  Compress  that  activity  until: 

(a)  either  Vcj  =  CUj  for  that  activity  (eg.:  &) 

(b)  or_  some  other  path  becomes  critical  (eg/  at  the 

i 

end  of  the  second  iteration,  ^4  is  reduced 
from  7  to  5  weeks  instead  of  3  weeks)  This 
situation  leads  to  multiple  critical  paths[eg/ 
(1,2,4)  and  (1,2, 3, 4)  for  a  project  duration 
of  9  weeks]. 
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( 1 v )  In  the  case  of  multiple  'vitical  paths,  step  (11) 
requires  the  computation  of  joint  cost  slopes,  [sg : 
at  the  end  of  the  second  Iteration,  the  joint  cost 
slope  of  (2,4)  and  (2,3)  Is  20  +  16=  36  and  that  of 
(2,4)  and  (3,4)  is  20  +  8  =  28].  Accordingly  choose 
minimum  among  the  set  of  cost  slopes  and  joint  cost 
slopes  relating  to  the  multiple  critical  paths  (eg 
choose  (2,4J  &(3,4) ) . 

(v)  Proceed  with  the  steps  (ii)  to  (iv)  until,  for  eacn 
activity  of  at  least  one  critical  path,  Y^j  36  GLij  and 
therefore,  no  more  compression  is  possible  [ecj[:  by 
compressing  both  (2,4)  and  (3,4)  simultaneously  by 

2  weeks  at  the  third  iteration,  the  activities  (1,2) 
(2,4)  are  now  set  at  their  crash  times  (e^  2  weeks  and 

3  weeks)  and  this  terminates  the  compression]. 

The  direct  project  cost  curve  for  this  example  is 

shown  in  Fig.  (4.11).  The  curve  is  piece-wise  linear. 

By  adding  to  this  an  indirect  cost  curve  (having  a  posi¬ 
tive  cost  slope)  a  convex  total  cost  curve  is  obtained. 


C-wr  PlggCT 


PROJECT  C05T  CUgV£ 


4.3  NETWORK  FLOW  APPROACH 


&y. 


Kelley  [65]  formulates  the  problem  given  above  as 
a  parametric  linear  programming  problem.  The  objective 
function  consists  a  maximization  of  the  project  utility 
which  is  a  function  of  the  project  duration.  The  con¬ 
straint  set  has  precedence  and  terminal  restrictions  of 
the  network  in  addition  to  the  upper  and  lower  bounds  on 
activity  duration.  Kelley  develops  a  fairly  efficient 
algorithm  by  Incorporating  a  slight  variation  of  Ford- 
F«jlkerson4s  Network  Flow  algorithm  [36]. 

In  a  companion  paper  by  Fulkerson  [39],  a  similar 
network  flow  approach  is  employed  to  minimize  the  total 
project  cost.  A  brief  summary  of  the  primal,  dual  and 
modified  dual  formulation  due  to  Fulkerson  is  given  in 
Fig.  (4.31).  As  illustrated  by  Elmaghraby  [121],  pp. 
(113-129),  Fulkerson's  flow  algorithm  may  be  divided  into 
five  steps.  Using  the  same  examr-1 2  given  in  the  last 
section,  details  of  repeated  computation  of  these  steps 
are  given  in  Fig.  (4.32).  Even  though,  the  modified  dual 

form  in  Fig.  (4.31)  indicates  the  need  for  two  parallel 
•  • 

arcs  for  each  (  C,j  ),  only  one  of  them  having  the  capacity 

C{,j  is  exhibited  in  Fig.  (4.32)  to  maintain  simplicity  of 

exposition.  The  second  arc  having  infinite  capacity  i r 

added  as  and  when  necessary. 

•  •  l 

Step  (1):  For  each  {  lj  J  )  and  (  K*  I  ,Z  ) ,  set  »  O  and 
Vi} ^  bij  i  compute  maximum  project  duration  (  A  ),  and 
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An  A')  ,(*} 

assign  label  {  Jij  }  Cij  -  TLj  )  or  (  fLi  f  oo  ) ,  the 
'sitter  one  for  the  Infinite  rapacity  arc. 

Step  (2):  Label  the  origin,  (j_e  node  1)  with  (c CjO). 

Check  each  arc  originating  from  node  1  for  the  relation- 
ship  Sjj  -  dij  -tj  ~0  ,  where  tj  is  the  earliest  oc¬ 

currence  of  event(  j  ) .  If  so,  label  node  (  j  )  with(oo^  I  ) . 
Proceed  to  label  successive  nodes  with  (oe^o),  if  the 
relationship  5ij  3  ti  +  Gtj  "  tj 3  0  holds  good.  If  the  ter¬ 
minal  node  (N)  can  be  labeled  (  «o, — ),  the  algorithm 
terminates.  Otherwise,  follow  the  subsequent  analysis, 
maintaining  the  node  labels  ( oO  ) . 

Step  (3);  Again,  start  with  the  origin  having  the  label 
(oOj  O  ) .  For  any  labeled  node  (  C  ) ,  check  al  1  nodes  (  j  ) 
connected  with  (c)  for  the  following  relationship 

Sij*  *  fcc+  blj  -  tj  -  O  a wL  <  Clj  •  5f  *o, 

/  .  r  /(O' 

label  node  (j  )  with  (  J  l  )  where  ^j=-mi  n  \  ("tj  J£j  ^ 
Should  the  relationship  j  *  ti  +  G£j -tj  =0  bold  good, 
label  node  (j)  with  ( A  )•  Continue  labeling .  If 
terminal  N  is  labeled,  a  break  through  has  occur  and 
so  proceed  to  step  (4).  Otherwise,  proceed  to  step  (5). 
Step  (4):  Start  with  node  (N),  For  any  labelled  node  (J  ) 
with  label  increase  the  flow  {  i'ij  )  by  (  ) » 

reduce  the  residual  capacity  in  arc  (  t,j  )  by  (  ^  )  and 

9 

erase  the  labeling  of  (J  ).  Continue  this  process,  until 
node  (1)  is  reached  and  go  back  to  step  (2). 


IMD  wmii)i)'ll  ||~|  nr  II 
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Step  (5);  Define  a  "cut  set"  Z=Cz»  j  Zzl  of  arcs  {<«>} 
such  that  one  node  Is  labelled  and  the  other  Is  not  labelled, 
where,  for  K*  i  , 

_(*>  ) 

Z{  a  (si)  |  (0  c*  labeled l  ,0)  U  not  UbeUA,  <°J 
Z  *  ^('si)lCO  ^  t^wtjJbeUA^O)  ^«-j 


■  l«0 


Ignore  all  arcs  of  Z|  where  the  equality  5j^  *0  holds. 

s,  =  rrg  Q.  Sif]  J  Scf  <  o 


Define 


Sa=  -gfrsn,  ^>o 

and  put  SsMin  <  S«  }  82, ) .  For  unlabelled  nodes,  change 
the  node  times  (  fcj  )  to  (  fcj  -  S  } .  Discard  all  node  labels 
other  than  {00 )  and  go  back  to  step  (2). 

This  iteration  would  ultimately  terminate  in  step  (2), 
when  the  terminal  node  (N)  is  labelled  (00)  indicating  that 
no  more  compression  of  project  duration  is  possible.  The 
direct  project  cost  curve  may  be  drawn  from  these  calculations. 


4.4  COMPUTATIONAL  SHORT  CUTS 

Some  original  simplifications  and  modifications  of 

Fulkerson's  algorithm  may  be  seen  in  an  article  by  Roper 

[102],  who  has  also  borrowed  ideas  from  Kelley.  Roper's 

algorithm  produces  sub-project  cost  curves  in  addition 

the  project  cost  curve.  Several  other  extensions  of  the 

network  flow  approach  of  Fulkerson  [39]  have  been  published 

by  different  authors.  The  list  includes  Pareikh  and 
* 

Jewell  [93]  who  provide  for  decomposition  of  large  net¬ 
works  into  severtisubnetworks  so  as  to  reduce  storage  demands 


*  Blanning  and  Rao  [93]  present  a  simpler  method 
which  does  not  require  the  use  of  Linear  Programming 
techniques . 
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in  machine  computation,  and  Handa  [50]  and  [51]  who  out¬ 
lines  computational  short  c*'ts  by  eliminating  network 
segments  which  do  not  Influence  the  project  cost  curves. 

In  a  different  content,  MacCrimmon  &  Ryavec  [79]  and 
Van  $lyke  [108]  suggest  similar  computatlonalshort  cuts 
(as  outlined  in  Section  3  of  this  paper),  while  Fulkerson 
[41]  &  [4-JL],  discusses  among  other  things,  an  iterative 
"longest  chain"  method  for  developing  the  project  cost 
curve  which  can  be  used  to  delineate  successive  critical 
paths. 

4.5  NODE  BALANCING 

If  the  convex,  continuous,  non-increasing  activity 
cost  curve  were  to  be  unbounded  at  both  ends  so  as  to 
resemble  a  hyperbola,  as  outlined  by  Clark  [18],  all 
activities  in  the  optimum  schedule  become  critical.  Berman 
[4]  relaxes  the  monotone  assumption  of  Clark  [18]  by  allow¬ 
ing  the  cost  curve  to  turn  upward  above  a  certain  activity 
duration.  Using  an  'evenion  node'  type  of  network  repre¬ 
sentation  where  several  activities  terminate  at  a  given 
node  while  several  others  start  from  the  same  node,  the 
nodes  are  balanced  one  by  one  by  equating  the  sum  of  time- 
cost  slope  of  incoming  activities  with  that  of  the  out¬ 
going  activities.  The  iterative  algorithm  terminates  when 
all  nodes  are  simultaneously  in  balance.  Ne  thus  derives 
a  minimum  cost  feasible  schedule  for  each  feasible  schedule 


for  each  feasible  project  duration  which  results  in  a 
continuous  convex  project  cost  curve.  He  applied  this 
principle  to  a  simple  serial  network  with  stochastic 
activity  duration.  He  observes  from  this  that,  In  gen¬ 
eral,  uncertainly  <n  activity  duration  tends  to  acceler¬ 
ate  the  occurrence  cf  intermediate  events.  In  particular, 
the  ones  which  precede  the  uncertain  activities. 


4.6  REPETITIVE  PROJECTS 

More  recently,  under  the  title  of  "Multicycle  Project 
Planning",  Fisher  and  Nemhauser  [35]  outline  a  method  of 
arriving  at  a  composite  project  cost  curve  Fig.  (4.63)  for 
a  serial  multicycle  project  Fig.  (4.61)  where  the  cost 
curve  for  each  cycle  Fig.  (4.62)  can  be  arrived  at  with 
the  well  established  methods  of  Kelley  [65]  and  Fulkerson 
[39].  In  real  situations,  the  repetitive  operations  are 
not  strictly  serial  since  portions  of  immediately  succeed¬ 
ing  cycles  can  begin  when  similar  activities  in  the  cur¬ 
rent  cycle  are  completed.  Such  semi-parallel  multi  cycle 
projects  lead  to  several  complications  in  applying  the 


time-cobt  trade  off  principles,  even  though  it  is  compara¬ 
tively  easy  to  identify  the  critical  path.  However,  they 
provide  an  algorithm  which  generates  the  jobs  and  prece¬ 
dence  relations  of  an  entire  multicycle  project.  If  a 
corresponding  description  of  a  single  cycle  is  given  in 
an  algebraic  form  (i£  a  set  of  Integers  representing  events 
and  a  set  of  ordered  pairs  of  integers  from  the  event  sets 
representing  the  activities). 


Fi 6.  (4>6£) 


Individual  CycU  Prolact  Cm*  Cwva, 


F~l  G  (|4 »  63^ 


MdUcydt  frajact  Cod  Cvrva. 
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4.7  STRUCTURAL  INTERPRETATION 

While  the  few  references  cited  above  relate  to  the 
extensions  of  Fulkerson's  Network  Flow  Algorithm,  for  the 
same  restrictive  assumptions  of  activity  cost  curve,  an 
entirely  different  method  of  computing  project  cost 
curve  is  given  by  Prager  [96].  The  use  of  Fulkerson's 
approach  necessitates  the  dual  form  of  the  orginal  Linear 
Programming  as  given  In  Fig.  (4.31),  but  Prager  works 
within  the  framework  of  the  original  problem  giving  a 
structural  Interpretation  to  flow  of  cost  in  the  network. 

4.8  USE  OF  INTEGER  PROGRAMMING 

The  restrictive  assumptions  of  continuity  and  con¬ 
vexity  inherent  in  Kelley,  Fulkerson  models  are  relaxed 
by  several  authors  [1],  [87],  and  most  of  them  do  not 
necessarily  yield  optimum  solutions.  However,  Meyer  and 
Shaffer  [85]  have  come  up  with  an  integer-linear  program¬ 
ming  technique  which  can  handle  such  cases  fc^*  networks 
having  fifty  or  less  activities.  As  pointed  out  by  Davis 
[24l  and  Carruthers  &  Battersby  [14],  the  integer  program¬ 
ming  approach  holds  great  promise  for  further  developments. 
In  the  areas  of  network  algebra,  cost  time  trade  offs  and 
resource  allocation,  recent  advances  by  Crowston  and 
Thompson  [22],  Balas  [2]  and  Pritsker  and  Watters  [99]  are 
based  upon  zero-one  integer  formulation.  More  references 
using  the  integer  programming  approach  are  cited  in  the 


sixth  section.  Crowston  has  further  extended  his  decision 
C.P.M.  concept  with  several  computational  short  cuts  to 
include  resource  allocation  models. 

Decision  C.P.M.,  (or  D. C.P.M.)  as  given  by  Croston  and 
Thompson  [22]  consists  of  a  formal  comparison  of  several  ad 

hoc  planning  networks  having  different  time,  cost,  resource 

den 

or  technological  depen^ies  so  as  to  choose  those  alterna¬ 
tives  that  minimize  the  total  project  cost.  The  auti ors 
recommend  the  repeated  use  of  the  same  method  during  the 
execution  phase  of  the  project. 

Using  activity-on-node  representation  of  a  network 
and  the  notations  given  in  Levy,  Thompson  and  Wiest,  Ch. 
(22),  [66],  they  define  a  decision  project  graph  G  as  a 
graph  with  nodes  representing  jobs  and  a  directed  time 
segment  connecting  two  nodes  Si,j  (in  job  set  5^)  and 
5mn(in  job  set  Sm ) ,  if  and  only  if  Sij  is  an  immediate 
predessor  of  Synr\  S<j«  5  wn.  They  define  a  de¬ 

cision  project  as  a  set  T  '  a  set  of  job  sets  S,  Szt  S*,- ••) 
together  with  the  specifi*-1  ..terdependencies  and  the 
relation  «  defined  on  o'.  The  outcome  of  decisions  based 
on  project  cost  analysis  would  be  to  eliminate  a  cluster  of 
jobs  from  the  decision  project  graph  to  obtain  the  final 
project  graph  used  in  a  regular  C.P.M.  approach.  This 
pruning  may  be  achieved  by  several  means  such  as  complete 
enumeration,  zero-one  Integer  Programming,  Heuristics  with 
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partial  enumeration  and  branch  and  bound  techniques. 
Crowston  and  Thompson  Illustrate  wi lh  a  few  examples 
the  several  ways  In  which  the  0-1  variable  can  be  used 
to  represent  complicated  Interdependencies  among  de¬ 
cisions.  The  exponential  growth  of  constraint  sets 
as  the  number  of  activities  are  increased  poses  several 
practical  problems  some  of  which  are  currently  being 
looked  Into  by  Crowston  and  others  at  M.I.T. 

4.9  C.P.M.  COST  VS.  C.P.M.  TIME 

While  the  papers  cited  above  under  "Cost-Time  Trade 
Off"  represent  the  mathematical  developments  in  this  area, 
the  doctoral  dissertation  of  Gessford  [43]  is  devoted  to 
a  comparative  study  of  CPM/cost  and  CPM/time.  Using  the 
data  derived  from  questionnaires  on  housing  construction 
as  well  as  simulation,  Gessford  concludes  that  the  utility 
of  CPM/cost  is  correlative  to  project  scope.  Hence,  med¬ 
ium  and  large  size  construction  firms  may  find  it  econom¬ 
ically  and  administratively  advantageous  to  add  cost  con¬ 
straints  to  their  existing  CPM/time  systems.  From  the 
point  of  view  of  Accounting  and  Finance,  CPM/cost  promotes 
cash  flow,  budgeting  and  analysis,  cost  estimation  and 
bidding,  cost  variance  control  and  analysis,  profit  plan¬ 
ning  and  the  integration  of  project  and  fiscal  planning. 

4.10  RELIABILITY  MATURITY  INDEX 

Finally,  Malcolm,  one  of  the  four  authors  of  the 
classic  paper  on  PERT  [80]  extends  the  basic  principles 


of  PERT  into  Reliability  Management  [81],  coining  a  new 
term  R.M.I.  which  represents  Reliability  Maturity  Index. 
This  descriptive  paper  outlines  the  procedure  to  be 
followed  in  measuring  reliability  in  addition  to  time 
and  cost  data  and  including  them  in  an  integrated  manage¬ 
ment  system.  The  system  makes  use  of  information  sampling 
of  management  reports  to  reduce  the  volume  of  input 
information. 


n 


SECTION  FIVE:  RESOURCE  ALLOCATION:  HEURISTIC  MODELS 

As  In  Section  (4),  the  review  articles  of  Rosenbloom 
[103l;Davis  [24]  and  Carruthers  and  Battersby  [14]  do 
good  justice  to  the  development  of  resource  allocation 
models.  Martino  [c3]  has  devoted  an  entire  book  to 
outline  various  heuristic  models  while  Moder  and  Philips 
[87]  present  a  few  representative  models.  In  all,  these 
books  and  review  articles  cover  most  of  the  work  done 
in  this  area  and  in  particular,  a  comparative  study  by 

Pascoe  [94]  and  by  Mize  [86]  of  different  resource  alloca- 

a 

tion  models  provides  a  sound  b^is  for  an  objective  analysis 
of  this  area  which  has  currently  become  a  focal  point  for 
further  research.  This  is  evidenced  by  recent  published 
articles  and  a  detailed  review  of  these  articles  is  given 
here. 

5.1  T?!E  COMBINATORIAL  PROBLEM 

The  time-cost  trade  off  principles  outlined  in  the 
last  section  assume  an  unlimited  availability  of  resources. 
Davis  [24]  recognizes  the  time-cost  trade  off  procedures 
as  one  of  the  three  special  cases  of  resource  allocation 
procedures,  the  other  two  being,  the  smoothing  of  resource 
usage  if  sufficient  resources  are  available  or  the  schedu¬ 
ling  cf  activities  when  the  resource  availability  is  re¬ 
stricted  to  a  predetermined  maximum.  Either  way,  as  Kelley 
[66]  points  out,  there  are  two  major  sources  of  difficulty-- 
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how  should  the  problem  of  resource  loading  be  formulated 
and  how  can  it  be  solved?  Irrespective  of  whether  the  re¬ 
source  required  Is  labor,  money,  equipment  or  facilities, 
the  existing  practices  in  Project  Management  Indicate  a 
lack  of  explicit  criteria  for  obtaining  optimal  use  of 
resources.  Hooper  [55]  regards  this  as  a  major  problem 
area. 

With  regard  to  the  mathematical  formulation  of  the 
problem,  difficulties  arise  because  of  its  combinatorial 
structure.  Kelley  remarks: 

"At  least,  the  several  models  we  have 
considered  are  completely  intractable.  However, 
the  difficulty  seems  to  be  intrinsic  to  the  prob¬ 
lem.  The  way  the  restrictions  on  the  sequence 
in  which  jobs  may  be  performed  interacts  with  re¬ 
source  requirements  and  availabilities  forces  a 
solution  set  which  is  ^connected.  This  property 
exists  even  if  we  assume  that  men  and  equipment 
are  infinitely  divisible. 

Problems  having  unconnected  solution  sets 
are  called  combinatorial.  Very  few  of  them  can 
be  solved  (in  the  mathematical  sense)  in  a  rea¬ 
sonable  time.  However,  this  does  not  imply  that 
useful  solutions  cannot  be  obtained.  Ad  hoc 
intuitive  methods,  especially  when  used  in 
conjunction  with  a  computer  to  do  the  leg  work. 
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can  be*  quite  effective  for  obtaining  satis¬ 
factory  solutions." 

Carruthers  and  rattersby  [14]  summarize  the  state 
of  art  as  follows: 

"No  algorithms  f>r  the  resource  smoothing 
problem  exist  which  are  rigorous,  general  and 
practicable;  some  exist  whicn  are  rigorous,  or 
rigorous  and  practicable  but  only  in  special 
cases.  The  rigorous  approach  to  resource  smooth¬ 
ing  considers  the  entire  project  as  one  entity, 
and  each  decision  in  scheduling  *s  made  against 
the  global  requirements  of  the  Droject.  However, 
all  practicable  algorithms  are  intended  to  be 
reasonable  rather  than  optimal,  with  regard  to 
an  objective  function;  consequently  they  are  based 
on  local  decisions.  Resource  smoothing  is  cur¬ 
rently  a  process  of 'me! iorization '  rather  than 
'optimization' ." 

5-2  CLASSIFICATION  OF  HEURISTICS 

Mostly  the  heuristic  models  can  be  classified  under 
series,  parallel  or  series  parallel  methods  of  scheduling 
resources,  as  suggested  by  Kelley  [66].  The  original 
series  method  consists  of  ordering  the  list  of  activi¬ 
ties  by  increasing  slack  within  precedence  constraints. 

By  scheduling  one  activity  at  a  time  in  the  serial  list, 
the  precedence  ordering  guarantees  that  each  task's  logical 
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predecessors  will  have  been  scheduled  previously.  Thus, 
if  sufficient  resources  are  available,  the  activity  In 
hand  is  scheduled  to  start  at  the  earliest  date,  which 
follows  the  latest  completion  of  the  predecessor  activity. 
In  the  parallel  method,  all  activities  that  are  ready  to 
be  scheduled  are  assembled  on  a  given  date.  Using  some 
intrinsic  or  explicit  consideration  as  the  case  may  be, 
each  activity  is  examined  and  the  best  among  them  are 
scheduled  if  resources  are  available.  The  ones  which 
are  not  scheduled  are  brought  forward  to  the  next  period 
and  examined  along  with  the  new  activities  which  become 
available  then.  This  operation  is  continued  until  all 
activities  are  scheduled. 

While,  in  the  series  method,  the  original  sequence  is 
maintained  throughout,  there  is  no  such  identification  of 
activity  listing  in  the  parallel  method.  The  series 
parallel  method  begins  by  scheduling  exactly  as  in  serial 
method,  but  the  sequence  of  activities  remaining  to  be 
scheduled  is  reviewed  at  each  step  and  altered  if  neces¬ 
sary,  recognizing  the  changing  circumstances  which  result 
from  scheduling  new  activities. 

5.3  A  RESOURCE  -  TIME  TRADE  OFF  MODEL 

One  of  the  early  articles  in  this  area  is  due  to 
McGee  and  Markarian  [84]  who  describe  a  serial  method 
with  variation  in  resource  assignment  and  provision  for 


rescheduling.  They  consider  a  multi  project  situation 
where  resources  are  allocated  to  different  projects  on 
the  basis  of  the  difference  between  the  expected  and 
desired  project  completion  dates.  They  employ  a  pro¬ 
cedure  very  similar  to  cost-time  trade  off  principle 
outlined  in  Section  (4).  The  desired  schedule  is  de¬ 
termined  by  an  iterative  technique  starting  with  minimal 
manpower  allocation  for  each  activity  and  successive 
manpower  increments  along  the  most  productive  critical 
activity.  Even  though  the  technique  is  attractive  in 
view  of  its  rather  unique  resource-time  trade  off  struc¬ 
ture,  the  lack  of  computer  programs  or  computational 
results,  the  approach  has  not  received  due  consideration 
by  other  researchers. 

5.4  SIMULATION  ALGORITHMS 

On  the  other  hand.  Levy,  Thompson  aid  W1est[74]  also 
consider  a  multi  project  situation  in  a  job  shop  with  an 
elaborate  computer  program  which  minimizes  the  peak  load 
over  a  future  time  interval,  treating  resource  upper 
limit  as  a  variable  rather  than  as  a  restraint.  A  novel 
feature  of  this  method  is  the  use  of  random  choice  steps- 
in  rescheduling  slack  activities,  thereby  necessitating 
repeated  trials  inherent  in  simulation  techniques,  so 
as  to  arrive  at  a  satisfactory  schedule.  As  Davis  [24] 
points  out,  this  does  not  always  produce  optimal  results. 


77  „ 

since  Its  effectiveness  will  vary  with  the  number  of  non 
critical  activities  and  their  slack  values. 

One  of  the  authors  of  the  article  mentioned  above, 

J.  0.  Wlest  [114],  takes  another  look  at  the  original 
problem  as  one  allocated  on  the  basis  of  a  group  of 
scheduling  rules  rather  than  as  a  resource  smoothing 
problem.  In  the  basic  model,  In  any  given  period,  the 
currently  available  activities  are  listed  according  to 
their  total  slackness  such  that  the  most  critical  activi¬ 
ties  have  the  highest  probability  of  being  scheduled 
first.  Using  Monte  Carlo,  as  before,  the  activities 
are  scheduled  up  to  the  resource  limit.  The  unscheduled 
jobs  are  pushed  forward  and  when  they  eventually  become 
critical,  they  receive  top  priority  in  the  current  ac¬ 
tivity  list. 

As  of  today,  Wiest's  model  is  the  most  promising  one 
In  view  of  Its  flexibility  to  Include  a  variety  of  realistic 
system  attributes.  For  instance.,  his  sub-routines  permit 
splitting  and  staggering  of  activities,  trade  offs  between 
time  and  resource  for  Individual  activities,  use  of  multi¬ 
ple  resources,  transfer  of  resources  between  activities  and 
pre-emption  of  activities  as  well.  The  model,  thus,  pro¬ 
ceeds  with  scheduling  and  rescheduling  each  day  optimizing 
locally  a  complicated  cost  function  based  on  resource 
utilization  and  expected  duration.  His  cost  function 
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Includes  costs  of  unused  resources  and  premium  and 
overtime  resources  which  are  affected  by  the  schedule, 
overhead  costs  related  to  project  duration,  hiring  and 
lay  off  costs  which  change  the  resource  levels  and  so 
on.  His  model  is,  in  part,  based  on  the  concepts  of 
"critical  sequence1  and  "conditional  slack'  derived 
from  a  limited-resource  modification  of  network  sched¬ 
uling  methodology,  as  given  in  [113]. 

5.5  OTHER  HEURISTIC  MODELS 

Being  the  most  recent  one,  Wiest's  model  represents 
a  substantial  improvement  over  several  heuristic  models 
in  resource  allocation.  Among  them  are  the  contributions 
by  Burgess  and  KilleLtew  [11]  who  use  the  sum  of  squares 
of  resource  requirements  as  a  measure  of  effectiveness 
while  DeWitte  [25]  minimizes  the  absolute  magnitude  of 
flunctuations  from  a  calculated  project  mean  level  of  re¬ 
source  usage. 

Other  approaches  at  multi  project  level  include  the 
work  of  Mize  [86]  who  developed  a  new  heuristic  algorithm 
for  scheduling  activities  to  several  operating  facilities. 
Martino  [83]  has  developed  a  heuristic  technique  called 
M.A  n.  to  consider  multi  project  resources  with  fixed  or 
variable  or  a  combination  of  fixed  and  variable  limits. 

His  technique  prescribes  no  upper  limit  either  to  the 
number  of  resource  types  or  to  the  number  of  activities. 
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Yet  another  system  available  in  the  market  is  due 
to  C.I.E.R.,  a  computer  service  organization,  who  offer  a 
comprehensive  package  caled  RAMPS  [71].  Being  a  proprietary 
system,  details  of  the  technique  are  not  available.  The 
system  is  believed  to  be  heuristic  and  performs  point-wise 
optimization  by  ranking  each  of  the  feasible  subsets  of 
currently  available  activities  on  the  basis  of  a  complex 
function  which  includes  work  continuity,  project  delay, 
idle  resource,  etc.  Another  recent  article  on  RAMPS  II 
by  Klein  [69]  provides  some  more  insight  into  the  critical 
Path  Algorithm,  describing  error  detection  and  network 
analysis  features.  It  can  operate  rapidly,  since  no 
topological  ordering  of  the  network  needs  to  be  performed. 

It  can  process  randomly  numbered  or  even  alphanumerically 
named  nodes. 

As  in  RAMPS,  several  computer  programs  Involving  the 
total  systems  approach  as  well  as  special  purpose  programs 
are  reported  at  the  29th  annual  meeting  ,0RSA,  Houston, 

Texas  (Nov.  1965).  Among  them  are  systems  such  as  d.P.S.S., 
due  to  Metzger  and  Thiele  and  SCOPE,  due  to  Miller  and 
Miller.  In  1964,  Moder  and  Philips  [87]  reported  the  de¬ 
tails  of  36  out  of  60  or  more  programs  then  available.  A 
similar  listing  of  currently  available  programs  would  in¬ 
deed  be  very  lengthy  and  perhaps  less  useful. 


5.6  COMPARATIVE  STUDIES 


Many  researchers  In  this  field  have  been  keenly  in¬ 
terested  in  a  comparative  study  of  different  procedures. 
For  instance.  Pascoe  [94]  has  examined  both  series  and 
parallel  algorithms  with  a  wide  range  of  characteristics 
involving  project  logic,  look  ahead  features,  time  and 
or  resource  constraints,  network  complexity  and  size, 
multiplicity  of  resources  and  with  five  different  objec¬ 
tive  functions.  His  conclusions  are  that  the  parallel 
methods  of  resource  smoothing  are  slightly  (ij;  at  0.1% 
probability  level  of  significance)  superior  to  the  serial 
method  and  tha!  the  use  of  the  latest  start  for  large 
networks  and  the  latest  finish  for  small  networks  are 
better  than  early  start  or  early  finish  routines.  Other¬ 
wise  there  is  hardly  any  difference  between  the  heuristics 
employed . 

Crowston  [23]  has  developed  several  ranking  criteria 
including  the  "greatest  remaining  resource  utilization", 
"remaining  slack",  etc.  and  suggests  several  ad  hoc  rules 
for  pre-empting  activities.  The  results  of  his  study  are 
still  inconclusive. 

Mize  [86]  has  also  made  a  comparative  study  of  se¬ 
lected  heuristic  algorithms.  The  number  of  test  cases 
constructed  by  him  were  too  small  to  demonstrate  con¬ 
vincingly  that  any  one  method  was  superior  to  another. 
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In  testing  and  evaluating  heuristic  models,  Wiest 
[114]  suggests  three  feasible  approaches:  laboratory 
tests  as  outlined  by  Pascoe,  field  tests  where  a  given 
model  is  applied  to  a  variety  of  real  world  scheduling 
problems  under  differeent  scheduling  environments  and 
constraints  and  f  1  nal  1  •*  comparative  tests  where  several 
computer  models  and  practical  methods  are  compared  for 
the  same  real  projects.  Much  work  needs  to  be  done  in 
each  case  before  objective  conclusions  can  be  reached. 

5.7  RESOURCE  SCHEDULING  FOR  STOCHASTIC  ACTIVITIES 

The  results  reported  so  far  refer  to  single  or  multiple 
projects  having  deterministic  activity  durations.  The 

e 

doctoral  dissertation  of  Fendl^y  [34]  on  the  other  hand 
concentrates  on  the  multi  project  resource  scheduling 
given  that  the  activities  have  stochastic  duration  of  a 
PERT  type.  Using  the  Monte  Carlo  Technique  on  a  combina¬ 
tion  of  hypthetical  projects  each  having  up  to  20  activi¬ 
ties,  three  resource  types  and  200  iterations  and  eight 
different  priority  rules,  Fendley  simulates  a  project 
scheduling  problem.  The  two  most  significant  conclusions 
from  the  experiment  are: 

1.  The  minimum-slack-first  (M.S.F.)  priority  rule, 
when  used  with  realistic  due  dates  is  the  best  of  those 
rules  tested. 

2.  Realistic  due  dates  may  be  set  by  an  analysis 
of  resource  load  on  the  facility  to  determine  the  amount 


of  slippage  that  must  occur  to  perform  all  projects 
with  fixed  resources. 

Under  the  M.S.F.  rule,  an  activity  is  loaded  from 
a  quelle  with  a  priority  dependent  upon  both  project 
characteristics  and  the  relative  urgency  of  the  project's 
due  date.  Using  regression  analysis  on  the  experimental 
data,  Fendly  predicts  slippage  necessary  for  the  per¬ 
formance  of  all  projects  in  the  hypothetical  system,  so 
that  the  system  can  be  loaded  to  any  desired  capacity. 

Even  though  it  is  hard  to  generalize  from  this  limited 
experimental  investigation,  the  approach  of  Fendley  is 
rather  unique  and  promising. 

5.8  PRACTICAL  STUDIES  IN  COMPUTATIONAL  FEASIBILITY 

Another  recent  doctoral  dissertation  by  Kapper  [122], 
who  applied  the  Critical  Path  Analysis  to  health  and 
hospital  fields,  provides  methods  for  handling  several 
practical  and  computational  difficulties.  Among  them 
are:  (a)  the  rapid  and  practical  segmentation  of  any  plan 
into  manageable  timesteps  (b)  the  comparison  of  the  re¬ 
source  scheduled  for  use  with  those  actually  available 
at  any  point  of  time  (c)  the  geryal  ization  of^schedule 
where  the  resources  which  are  used  approach  but  not  exceed 
the  resource  limits  specified  for  each  point  of  time  (d) 
the  rapid  practical  integration  of  this  algorithm  and  its 
associated  computer  programs  with  other  currently  available 


scheduling  programs  (e)  a  variety  of  distribution  functions 
for  characterizing  activity  completion  times  (f)  the  bias 
associated  with  time/cost/resource  requirement  estimates 
(g)  the  production  of  an  on-line  real  time,  planning- 
scheduling-all  ocation  capabili ty  consistent  with  the  con¬ 
tingency  and  feed  back  requirements  of  plan  implementa¬ 
tion;  and  (h)  the  incorporation  of  auxiliary  capabilities 
(e.g.  decision  table  logic,  trend  analysis  and  information 
retrieval)  in  order  to  achieve  sub  optimal  allocation 
(Thesi s  Abstract) . 


SECTION  SIX:  RESOURCE  ALLOCATION 


The  papers  discussed  so  far  in  the  area  of  resource 
allocation  have  been  mostly  heuristic,  with  particular 
reference  to  elaborate  schemes  outlined  by  Wiest  [114], 
Fendley  [34]  and  Kapper  [122].  However  the  most  recent 
work  in  this  area  is  directed  towards  a  mathematical 
formulation  of  the  problem.  In  brief,  the  review  in¬ 
cludes  a  significant  contribution  by  Johnson  [63]  who 
applies  the  principles  of  branch  and  bound  for  simple 
models,  the  work  of  Pritsker  and  Watters  [99]  who  sug¬ 
gest  a  zero-one  programming  approach,  a  resource-time 
trade  off  model  of  Ghare  [44],  an  excellent  series  of 
Russian  and  Yugoslav  papers  [109],  [12],  [95]  covering 
the  areas  of  quadratic  programming ,  dynami c  programming 
and  the  optimal  control  theory.  This  section  also  con¬ 
tains  the  dissertation  abstract  of  Bennett  [123]  who 
reports  an  Integer  Programming  approach  and  a  series  of 
con ‘■ri jutions  relating  the  assembly  1  ine  balancing  to  the 
resource  allocation  problem. 

6.1  BRANCH  AND  BOUND  IN  RESOURCE  ALLOCATION 

The  branch  and  bound  idea  was  used  by  with  Little, 
et  al  [75]  while  they  were  seeking  a  manageable  solution 
for  the  travelling  salesmen  problem.  Lawler  and  Wood  [73] 
provide  a  survey  of  several  papers  in  this  area,  illustrat¬ 
ing  the  wide  applicability  of  this  concept  in  several  O.R. 


problems.  Johnson's  thesis  [61]  is  yet  another  illustra¬ 
tion  of  this  simple  but  powerful  method  of  branch  and 
bound.  Using  both  critical  path  and  resource  require¬ 
ment  data  he  outlines  a  minimum  bounding  technique  which 
automatically  discards  several  partial  schedules  which 
cannot  possibly  improve  upon  a  known  complete  schedule. 

He  defines  such  an  elimination  procedure  In  terms  of  a 
“limited  enumeration"  algorithm.  He  also  provides  for  a 
comparison  among  partial  schedules  themselves  thereby 
reducing  substantially  more  computational  work 

These  "partial  schedule  dominance"  tests  are  the 
added  features  of  his  computer  program  which  is  workable 
up  to  50  jobs.  Beyond  this  range,  the  computational  time 
required  become  too  large  to  be  economical.  However,  the 
branch  and  bound  approach  may  still  provide  a  meaningful 
approximations  to  large  complex  problems.  Depending  upon 
the  nature  of  the  problem,  the  use  of  appropriate  bound¬ 
ing  and  dominance  tests  along  with  an  efficient  data 
organization  should  considerably  improve  upon  the  pre¬ 
sent  state  of  the  art. 

6.2  A  ZERO  -  ONE  INTEGER  PROGRAMMING  FORMULATION 

Pritsker  &  Watters  [V3]  present  a mathemati cal  approach 
to  scheduling  problems  in  general  with  the  help  of  a  zero- 
one  integer  programming  formulation.  They  demonstrate  that 
it  can  accomodate  a  wide  range  of  real  world  situations 


and  perhaps  large  practical  problems  may  eventually  be 
solved  in  this  manner.  In  particular,  the  0-1  variables 
are  used  to  indicate  for  select  periods  whether  or  not 
a  job  has  been  completed.  The  period  selection  depends 
upon  the  job  arrival  time,  due  dates,  sequencing  relation¬ 
ships,  etc.  They  consider  three  objectives:  minimize 
total  throughput  time  for  all  projects;  minimize  the  time 
by  which  all  projects  are  completed  and  finally  minimize 
total  lateness  or  lateness  penalty  for  all  projects.  The 
constrained  set  includes  limited  resources,  precedence 
relations  between  jobs,  job  splitting  possibilities,  pro¬ 
ject  and  job  due  date«,  resource  substitution,  concurrent 
and  non  concurrent  job  performance  requirements.  With 
the  help  of  simple  theoretical  examples,  they  demonstrate 
the  superiority  of  0-1  formulation  over  known  heuristics. 
This  attractive  theoretical  formulation  remains  to  be 
tested  for  large  practical  problems  with  regard  to  its 
computational  feasibility  even  as  Johnson's  "branch  and 
bound"  [63],  Meyer  anJ  Shaffer's  Cost-Time  Trade  Off  [85] 
and  Crowston  and  Thompson's  Decision  C.P.M.  [22]  have  all 
encountered  similar  difficulties.  Yet  another  theoretical 
approach  to  resource  allocation  problem  using  integer¬ 
programming  may  be  seen  in  Hadley  [49]. 

6.3  ANOTHER  RESOURCE  -  TIME  TRADE  OFF  MODEL 


As  suggested  by  Clark  [18],  Berman  [4],  Jewell  [60] 
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and  McGee  and  Markarian  [84],  the  work  of  Ghare  [44]  is 
based  on  the  assumption  that  the  duration  of  an  activity 
can  be  controlled  by  the  modification  of  resources  al¬ 
located  to  the  activity.  Time  and  Resource  are  con¬ 
sidered  the  through  variable  and  the  across  variable  of 
the  allocation  system  respectively.  Using  the  techniques 
of  the  systems  theory,  and  a  mathematical  programming 
algorithm,  the  project  duration  is  minimized.  The 
major  emphasis  Is  on  the  management  control  of  activity 
duration,  while  proceeding  with  an  optimal  allocation  of 
scarce  resources. 

6.4  A  QUADRATIC  PROGRAMMING  APPROACH 

Significant  among  the  analytical  contributions  on 
Resource  Allocation  are  the  works  of  Voronov  and  Petrushinin 
[109]  who  use  Beale's  Quadratic  Programming  method  in  mini¬ 
mizing  simultaneously  the  mean  square  deviation  of  the  re¬ 
source  allocation  function  from  a  constant.  The  allocation 
of  resources  Yij  (t)  to  activities  may  become  non- 

uniform,  changing  by  jumps  at  certain  time  intervals  while 
remaining  constant  in  between,  satisfying  the  condition 

f  sO  /  -fay 

Ycj  °)  -  <  >0  .  I  ^  i 

l  =°  j  o<  >  J 

where  fcerf  is  the  time  of  occurrence  of  event  o< . 

The  model  includes  the  possibility  that  the  allocation 


88. 


function  ^tjvv  reduces  to  zero  during  any  of  the 
intermediate  intervals.  An  important  assumption  of 
the  model  is  that  the  moments  at  which  events  occur 
are  assumed  to  be  fixed  and  consequently,  the  interval 
between  event  times  (  ~  )  are  also  fixed.  There 

fore,  for  every  activity,  there  is  a  linear  algebraic 


equation 

£  Vq(^+0)  =  f„ 

o^*0 


(t^j)  and  (  i/j  -  l  ,1, . T) ) ,  where  QCj  is  the  initially 

given  work  content  of  the  activity.  The  objective  function 
has  a  quadratic  form 

IT  -  £  Vi.j C 


ot-o  — 1 

where  the  activity  is  a  member  of  the  subset  of 

the  work  performed  by  resources  of  the  S**  type.  With 
the  help  of  a  simple  theoretical  example,  Voronov  and 
Petrushinin  show  all  the  steps  involved  in  using  Beale's 
Method  of  Quadratic  Programming. 

6.5  USE  OF  OPTIMAL  CONTROL  THEORY 

Burkov  [12]  considers  a  class  of  resource  distribu¬ 


tion  problems  from  the  standpoint  of  the  optimal  control 
theory.  He  regards  the  process  of  accomplishing  a  pro¬ 
ject  as  a  transient  process  in  the  automatic  control 


I 

I 


system  where  the  resources  play  the  role  of  control 
actions.  The  dimension  of  phase  space  of  the  system  state 
is  equal  to  the  minimal  number  of  independent  paths  of  the 
project  network  (i_e  the  paths  having  no  common  nodes). 
Burkov  confines  his  attention  to  the  problem  of  cost 
minimization  within  a  restricted  project  completion  time. 
He  shows  that  the  convergence  problem  in  the  iteration 
process  involved  in  searching  for  an  optimal  solution 
is  reduced  to  analysis  of  system  stability  in  a  transient 


process  of  automatic  control  system. 

6.6  A  DYNAMIC  PROGRAMMING  FORMULATION 

Burkov's  control -theoretic  approach  has  been  limited 
to  cost  minimization.  Pursuing  his  ideas  further,  Petrovic 


[95]  maps  the  predecessor-successor  relations  of  activities 

€  A  [  into  a  set  of  transforms  , 

l  1  ,  i  ns  0,1,.  .-  N-l 

Xcj  O+0  ■  T<J  \  M  ,  Uy  (n)  |  ,  €  A 

where  the  state  variable  is  the  scope  of  work  Xcj  C 0 
required  to  complete  an  activity  and  the  decision  variable 
is  the  resource  (JLcj  (t) allocated  to  each  activity.  If 
the  availability  of  a  resource  is  limited  to  U,  then 
,2L  ,  "h -  0,ly- -*N-I.  The  admissible  set  U  may 

be  regarded  as  a  given  function  of  time  or  the  state  of 


the  project  itself.  Further  using  the  network  constraints 


Uj<<. 

& 

*9  *a  1 

1  (cjj)  6  * 

v  r3'  A 

Ujfc  -0 

ftrf 

Xj+c  =o  . 

Cj,<0  6  A 

9  *  0 

CVnoi  £  ^  j  ^ 

and  scalar  objective  function, 

t=  z.  Y_  C  (t');  Uijw) 

Me  a  »-o 

Petrovic  views  this  multistage  decision  process  as 
a  dynamic  programming  problem  having  the  following  re¬ 
currence  relationship  in  the  case  of  a  single  resource. 


■fN.K  (xyt*-*-  «-*♦&) 


Mid 


-  r,  N  ,  -n  \  L  J  J  l. 

Uij  -*<),  (N-K+yj  ^ 

frr  k «  2,  3, -  N  ^ 

wlth  f  /  o 

it-l  [X'-j  V'^j  Ui:fe(N-l,N)L  -J 

Inevitably,  the  dimensions  of  the  functional  equations 
which  appear  in  the  model  become  very  high  in  practice, 
thereby  presenting  serious  limitations  in  computational 
feasibility.  In  reducing  the  necessary  amount  of  computa¬ 
tional  effort,  Petrovic  suggests  three  sophisticated 
approaches.  The  first  approach  is  called  sub-project 
programming  which  involves  an  ingenious  time-wise  decompo¬ 
sition  of  a  time-scale  network.  Selecting  sub-projects 
in  this  manner,  every  crossing  state  variableixjtween  two 
adjacent  sub-projects  is  considered  as  an  output  of  one 
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of  them  and  an  input  of  the  other.  An  optimum  aggregation 
of  all  sub-projects  can  then  be  made  vnth  the  help  of  the 
crossing  variables  in  such  a  way  that  the  aggregate  solu¬ 
tion  is  equivalent  to  that  of  the  original  network. 

The  second  approach  is  again  another  form  of  decom¬ 
position  employing  a  'two  level  concept'  which  is  based 
on  the  results  of  Lasdor.  and  Schoefflar  [72].  The  two 
level  scheme  consists  of  the  decomposed  component  sub¬ 
systems  to  each  of  which  is  assigned  a  specially  con¬ 
structed,  independent  sub-problem.  The  objective  func¬ 
tion  of  each  sub-problem  contains  a  number  of  variable 
parameters  or  "prices"  whose  values  are  determined  by  a 
second  level  controller.  Sub-problem  solutions  are 
functions  of  these  parameters.  Lasdon  and  Schoeffler 
give  the  criteria  which  enable  the  second  level  to  find 
the  "optimal"  parameter  set  by  a  simple  iterative  procedure. 

The  last  approach  represents  an  approximation  in 
policy  space  as  given  in  Bellman  and  Dreyfus  [3].  The 
procedure  consists  of  a  simple  initial  approximation  to 
the  optimal  policy.  Then  a  control  policy  satisfying 
certain  properties  is  obtained  such  that  it  yields  a 
better  performance.  The  control  policy  is  then  changed 
partially  to  effect  optimal  improvement  at  each  iteration. 
When  the  value  of  the  criterian  decrement  becomes  less 
than  the  specified  tolerance,  the  iterations  may  stop. 
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even  though  computations  may  be  stopped  much  earlier  if 
desi red . 

Petrovic's  paper  represents  a  significant  departure 
from  the  familiar  heuristic  approach  to  the  resource  alio 
cation  problem.  However  the  ideas  presented  here  are  yet 
to  be  validated  by  specific  illustrations  which  would 
demonstrate  their  computational  feasibility. 

6.7  COMPUTATIONAL  EXPERIENCE  WITH  DIFFERENT  APPROACHES 

Yet  another  integer  programming  approach  to  the  re¬ 
source  allocation  problem  is  given  in  an  unpublished 
doctoral  thesis  of  Bennett  from  Cornell  University  [123]. 
An  excerpt  from  his  thesis  abstract  is  given  below. 

"This  thesis  develops  several  approaches  to  the 
following  problem:  Sechdule  the  activities  of  a  pro¬ 
ject  network,  some  or  all  of  which  require  resources 
which  are  available  in  limited  amounts,  in  such  a 
way  that  the  activities  are  performed  in  proper 
sequence,  the  resource  availabilities  are  not  ex¬ 
ceeded,  and  the  duration  of  the  project  is  as  short 
as  possible. 

Mathematical  models  are  presented  which,  when 
solved,  can  assure  that  these  conditions  are  met. 

A  linear  model  represents  the  case  in  which  activity 
durations  are  considered  to  be  fixed,  and  the  re¬ 
sources  of  constant  amounts  are  assigned  to  each 


activity.  A  linear  model,  a  quadratic  model,  and 
a  transportation  format  ranking  procedure  are 
given  for  the  situation  in  which  activity  durations 
may  vary  for  different  resource  assignments  and 
only  total  resource  requirements  are  given  for 
each  activity. 

The  two  linear  models  require  integer  linear 
programming  for  solution.  Experience  in  attempting 
to  solve  the  formulations  with  two  different  com¬ 
puter  codes  on  a  Control  Data  1604  computer  is 
reported.  Although  this  experience  was  valuable 
in  terms  of  testing  a  rather  new  mathematical  tech¬ 
nique,  it  is  clear  that  this- approach  could  not  be 
applied  to  the  planning  of  actual  projects  because 
of  limited  computer  memory  capacity  and  excessive 
run  time.  For  reasons  that  are  not  clearly  under¬ 
stood,  all  attempts  to  solve  the  quadratic  model  by 
a  method  suggested  by  Wolfe  were  unsuccessful.  The 
final  technique,  a  ranking  procedure,  is  limited  to 
only  narrow  applications. 

Because  of  the  limitations  of  the  minimum  time 
models,  an  heuristic  algorithm  is  given  for  the  case 
in  which  activity  durations  may  vary.  This  program 
essentially  works  through  a  list  of  activities,  as¬ 
signing  resources  to  them  for  each  time  unit  of  the 


project  until  all  resources  aviilable  during  that 
time  unit'  are  assigned  or  until  no  more  activities 
can  be  scheduled.  Consideration  then  turns  to  the 
next  time  unit,  and  the  procedure  is  repeated. 

When  all  activities  have  been  assigned  their  total 
resource  requirements,  the  project  is  considered 
complete,  and  the  program  terminates.  The  result¬ 
ing  schedule  is  feasible,  since  resource  avail¬ 
abilities  are  not  exceeded,  but  not  necessarily 
optimal,  since  the  duration  may  not  be  as  short 
as  possible. 

The  program  can  schedule  projects  containing 
as  many  as  700  activities  and  having  durations  as 
great  as  400  days.  Seven  resources  may  be  scheduled 
simultaneously.  In  addition  to  the  total  require¬ 
ment  for  each  resource,  the  input  for  each  activity 
includes  the  maximum  amounts  of  each  resource  that 
can  be  assigned  during  any  time  period  and  the  mini¬ 
mum  amount  that  can  be  assigned  if  the  amount  is 
greater  than  zero.  Experience  with  the  program 
indicates  that  it  can  be  applied  successfully  in 
scheduling  reasonably  large  projects." 
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6.8  ASSEMBLY  LINE  BALANCING  PROBLEMS 

Brief  descriptions  of  resource  allocation  models 
of  Wilson  [115]  and  Black  [7]  may  be  seen  in  the  review 
article  by  Davis  [24].  Apart  from  suggesting  a  dynamic 
programming  approach,  Wilson,  for  the  first  time,  provides 
a  linkage  between  resource  allocation  problems  and  as¬ 
sembly  line  balancing  problems  as  outlined  by  Hu  [59]. 
Proceeding  on  similar  lines.  Black  [7]  adapts  Gjtjair- 
Nemhauser  algorithm  for  assembly  line  balancing  [48]. 

More  recently,  the  similarities  between  these  two  tech¬ 
niques  are  established  by  Moodie  and  Mandeville  [88]. 

Their  main  source  is  Bowman's  [9]  presentation  of  an 
Integer  Programming  method  of  assembly  line  balancing. 
Moodie  and  Mandeville  indicate  the  basic  differences  be¬ 
tween  the  two  problems  and  illustrate  their  approach  by 
obtaining  a  computer  solution  to  a  simple  theoretical 
problem.  The  merits  of  heuristics  for  multiproject 
resource  balancing  are  also  discussed. 


SECTION  SEVEN:  DIRECTIONS  OF  RESEARCH 


In  one  of  his  early  papers  (1961)  describing  the 
primal-dual  algorithm  for  the  time-cost  trade  off  prob¬ 
lem,  Kelley  [65],  p.  314,  remarked: 

"Quite  often  large  portions  of  a  project 
diagram  form  individual  projects  in  their  own 
right.  They  connect  into  the  whole  project  at 
only  two  points,  their  origin  and  terminus,  re¬ 
spectively.  Each  of  these  sub-projects  can  be 
solved  separately,  a  spectrum  of  solutions  and 
utility  functions  being  obtained  in  each  case. 

Each  sub-project  may  now  be  replaced  in  the  pro¬ 
ject  by  a  single  activity  with  the  same  utility 
function  as  the  sub-project  it  replaces.  The 
original  project  is  now  smaller  by  many  activi¬ 
ties  and  events  and  can  be  solved  more  easily 
than  before." 

The  work  of  Parikh  and  Jewell  [93]  (1965)  on  decom¬ 
position  of  project  networks  essentially  achieves  this 
goal,  using  Fulkerson  s  [39]  algorithm.  While  providing 
a  solution  to  a  project  consisting  of  sub-projects  which 
are  connected  in  series  or  parallel  or  a  combination  of 
both,  they  recognize  that  the  "generalized  clustered 
network"  is  more  difficult  and  no  simple  algorithm  is 
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available.  They  suggest  the  use  of  "out-of-kilter" 
method  [37]  for  each  Interlink  and  each  desired  project 
duration.  In  >  different  context,  Jewell  [61]  has  made 
use  of  network  flow  approach  to  Include  "divisibility" 
of  activities,  suggesting  the  need  for  techniques  in 
the  area  of  mul ti -commodity  network  flows  if  two  or  more 
activities  are  divisible.  He  has  also  hinted  the  pos¬ 
sibilities  of  extending  this  model  to  time-cost  trade 
off  area.  More  recently,  Jewell  [62]  has  made  further 
contributions  to  the  "out  of  kilter"  algorithm  while 
progress  is  being  made  in  the  area  of  multi-commodity 
network  flow  [125]  and  project  decomposition  [78]  ex¬ 
tending  the  computer  program  to  PERT  type  of  networks. 

Yet  another  stream  of  thoughts  pertaining  to  stochastic 

to 

activity  duration  relates^the  chance  constrained  program 
[16]  and  directed  sub-dual  algorithm  [15]  for  time  calcu¬ 
lations  only,  leaving  room  for  dev*lor'ment  in  the  time- 
cost  trade  off  area.  The  synthesis  of  these  ideas  may 
also  include  the  feedback  stochastic  networks  of  Elmaghraby 
[32],  [33],  Pri trker  and  others  [97],  [98],  [112],  [117]. 
Finally  the  developments  in  "Decision  C.P.M."  due  to  Crowston 
and  Thompson  [22]  could  be  merged  with  others  providing  for 
extensions  in  decomposition  of  integer  programs  [124], 
branch  and  bound  applications  [63],  the  use  of  filter  method 
[2]  and  disjunctive  constrains  [2]  and  [104]. 
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Kllngel  .  [70]  observer,  that  while  resource  alloca¬ 
tion  programs  reduce  the  deterministic  completion  time 
by  increasing  the  parallelism  of  the  network,  the  prob¬ 
ability  of  project  completion  by  the  new  date  is  also 
reduced.  On  the  other  hand,  Wiest  [113]  shows  that  the 
usual  concepts  of  'critical  path'  and  'job  slack'  lose 
their  meaning  in  the  scheduling  for  large  projects  with 
limited  resources.  He  suggest  the  use  of  a  concept 
called  'The  Critical  Sequence'  which  identifies  a  sequence 
of  critical  jobs  composed  of  either  a  technological  se¬ 
quence,  or  a  resource  sequence  or  both.  In  a  similar 
manner.  Van  Slyke  [108]  introduces  the  concept  of  "criticality 
index'  while  Noettl  and  Brumbaugh  [91]  introduce  information 
theoretic  concepts  in  network  planning.  In  the  light  of 
these  new  concepts,  and  the  mathematical  developments  using 
control  theory  [12],  quadratric  programming  [109]  and  integer 
programming  [99],  [49],  [123]  for  resource  allocation  prob¬ 
lems,  further  strides  cau  be  made  to  integrate  time-cost 
trade  off  with  resource  limitations  and  perhaps  a  new  di¬ 
mension  involving  performance  standards.  Equal ly  .chal¬ 
lenging  are  the  fast  developing  areas  of  assembly  line 
balancing  and  production  scheduling  where  the  similari¬ 
ties  with  the  project  networks  have  already  been  estab¬ 
lished:  [115],  [7],  [88]  and  [107]. 


From  a  practical  point  of  view,  the  contributions 
of  King  and  others  [67]  and  [68]  need  to  be  substantiated 
by  further  data  about  the  estimating  behavior  of  indivi¬ 
duals  while  a  meaningful  synthesis  of  Line  of  Balance 
and  PERT  networks  (PERT/LOB)  as  suggested  by  Schoderbek 
and  Digman  [105]  extend  these  concepts  to  repetitive  jobs. 

In  conclusion,  the  exhibit  IV  attempts  to  provide  a 
basis  for  cross  fertilization  among  the  four  areas  of 
network  algebra,  P'iRT  assumptions,  cost-time  trade  off 
and  resource  allocation.  The  applicability  of  different 
mathematical  programming  and  other  techniques  in  this 
area  indicates  that  developments  in  these  respective 
areas  may  have  an  influence  upon  the  project  networks 
and  their  computer  programs. 
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u.  abstract  paper  brings  together  different  aspects  of  the  Critical  Path 

Analysis  exclusively  in  terms  of  the  mathematical  developments  in  this  area 
during  the  last  eight  years.  The  survey  is  divided  into  seven  sections. 

Tne  first  one  deals  with  the  four  stages  of  development,  each  successive 
stage  being  more  representative  of  the  real  situation.  The  next  five 
sections  give  explicit  reference  to  the  different  articles  which  contri¬ 
bute  to  the  stages  of  development  described  in  the  first  section.  The 
articles  are  reviewed  in  terms  of  the  network  algebra,  PERT  assumptions, 
cost/time  trade  off,  resource  allocation  (heuristic)  and  resource  alloca¬ 
tion  (analytical) ,  Finally,  the  seventh  and  last  section  of  this  paper 
provides  a  summary  of  ideas  and  possible  directions  of  research,  some  of 
•which  are  currently  being  undertaken. 
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